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ABSTRACT 


The use of central differences on a rectangular net, in solving the primitive or vorticity equations, produces 
solutions on each of two lattices. By exploring this lattice structure, a formal equivalence is established between 
the central-difference vorticity and primitive equations. A demonstration is given also that exponential instability 
previously found to result from certain types of boundary conditions is suppressed by applying these conditions in 
such a way as to avoid coupling the lattices. 


1. INTRODUCTION Partition the axes uniformly as follows: t/At=0, 1,2,..., 
and z/Ar=0, 1, 2,. . . , p, so that L=pAz. Using cen- 


Ina recent study (Platzman [7]) the primitive equations tral Gilérenves, 


(that is, the momentum and continuity equations) were 


integrated numerically in an investigation of surges on shou Pte, =) 
lake Michigan. It was pointed out there that, when where ¢=cAt/Az and 6 denotes the operator ° 

central differences are used, the finite-difference primitive 

equations give rise to solutions on two independent 5 u(x,t) =u(z,t+ —u(z,t— At) 


lattices. The purpose of the present note is to demon- 
strate that a more satisfactory understanding of the $A2,t)—u(2— 

analytical properties of these equations is provided by Subject to boundary conditions u(0,t)=0 and u(L,t)=0, 
taking account of this lattice structure.* we may represent u(z,t) by the discrete Fourier series 


2. THE WAVE EQUATION -1 ‘ 
Soult) sin kw (2.2) 
Consider a canal of uniform depth D and uniform cross 4 
section, closed by rigid vertical walls at z=0 and z=L. l=1. 2 at 
The differential equation for “tidal” disturbances is the 
one-dimensional wave equation where k, is the wave number of the /th mode. 
Substitution of (2.2) in (2.1) yields the following differ- 
ence equation for the expansion coefficient w,(t): 


—4 (o sin $k,Az)*u,. 
where is the particle velocity along z, and c=ygD. (o sin 


{This investigation was supported by funds provided by the U. 8. Weather Buresu. The general solution of this is 
Pa the substance of the present work was completed, the writer’s attention was 
‘sled to an elegant study by A. Eliassen, in which similar ideas are involved (A. Eliassen, =a cos v sin py 2.3 
A Procedure for Numerical Integration of the Primitive Equations of the T'wo-Parame- ult) - t+B t, ( ) 
‘er Model of the Atmosphere,” University of California at Los Angeles, Department of ys A 
Meteorology, 1956). where a and # are constants of integration, and 
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sin $v,At=o sin $k,Az (2.4) We first represent u(z,t) in the manner of (2.2), anj t 


determines the frequency. h(a,t) by the discrete series 


For initial conditions we specify wu and du/dt; in the sal 
finite-difference frame, this is equivalent to a specification h(z,t) = cos k yz. (3.2) 
of u,(0) and u,(At), in terms of which (2.3) becomes 
u,(t)=(ese v,At)[—u,(0) sin »,(t—At)+wu,(At) sin vf]. Substitution in (3.1) yields 
If we assign to u,(At) the value : 
6u,= +2(¢ sin $k Az)h, 


u,(At)=u,(0) +%,(0) At, 5,h;=—2(o sin $k,Az)u, 

where %,(0) is Ou,/ot at t=0, then we get for which we have the general solution 
u,(t)=(sec }v,At)[u,(0) cos v,(t—4At)+ u(t)=a cos vt+ sin 
1,(0)(¥At esc $v,At) sin vf]. h,(t)=B cos vt—a sin vt, 


with vy, as in (2.4). If, in addition to the initial distriby. 
tion of u, one specifies the initial distribution of h (instead 
of the initial acceleration 0u/dt), the solution is 


This form shows clearly the amplitude and phase distor- 
tion produced by truncation errors (see Platzman [5]), and 
expresses the solution explicitly in terms of data prescribed 


along the initial line. u,(t) = (sec 3v,At)[u,(0) cos +h, At) sin (3.38) 
The finite-difference system has 2(p—1) degrees of free- (¢)— (sec 4y,At)[h,(#At) cos vt—u,(0) sin 
dom, represented by the p—1 independent values of (3.3b) 


u(z,0) and the p—1 independent values of u(z,At), or by 
the values of the corresponding Fourier coefficients. 
These degrees of freedom are provided by the 2(p—1) 
distinct fundamental modes exp (+%»;t) determined by 
the frequency equation (2.4), of which p—1 correspond 

- to wave propagation along positive z, and p—1 correspond 
to backward propagation. 


Note that the initial values of A are those assigned on the 

line t=4At, because in the lattice of figure la, values of t 

h are not given on t=0. d 
In figure 1b is shown the lattice which is initiated by 

assigning h on the line t=0 and u on t=4At. The solution 

in this case is 

ur(t)= (sec cos sin 


3. THE PRIMITIVE EQUATIONS (3.4a) ‘ 
The preceding problem may be formulated in terms of hi(t)=(see $v,At)[hi(0) cos v4(t—$At)—w,(ZAt) sin vf) 
the linearized primitive equations (3.46) 
The lattices of figures la and 1b differ only by a phase Bs 
ou oh shift of }A¢ 
ot ox The equations (3.1) involve p—1 values of u (exclusive Bo 


of the boundary values, which are fixed at zero), andp § d 


oh =—D values of h; hence, there are 2p—1 degrees of freedom 

in the initial conditions. In order to obtain an equal 
where A is the free-surface displacement from mean level. number of fundamental modes for (3.1), one must include By 

As before, we place wu at the points t/}At=0, 2,4,..., the root 1=0, »,=0 of the frequency equation (2.4); the 

and 2/}Az=0, 2, 4,..., 2p; further, we place h at the corresponding mode is u=0, h=ho=constant, which ev- 


“half-way” points t/}At=1, 3, 5,... , and 2/$Az=1, 3, dently satisfies (3.1) and boundary conditions. Note d 
5, ..., 2p—1 (see fig. 1a). Then the central difference also that if (3.1b) is summed over all values of x where h iY 
is defined, the right side sums to zero, with the help a Bg, 


equations are, in the notation of (2.1), ae 
the boundary conditions, so that the z-sum of h is inde 


5u=—ab,h (3.la) pendent of t. (This is merely a statement of the integral 
i 5 h=—06,u; (3.1b) form of the continuity equation.) From (3.2), we sé 
that the z-sum of h is equal to pho; hence, again, io=co0- 

in (3.1), wand h have been reduced to dimensionless form stant. Sj 


through division, respectively, by ygD and D. If one To insure that the finite-difference solution is stable, it By, 
eliminates h from (3.1), by applying 6, to (3.1a), equation is sufficient to choose At and Az so that o =cAt/Ar<l, fir 


(2.1) is obtained, so that our frequency equation is again since then the frequency »; determined by (2.4) must be bs 
(2.4). However, instead of working with the wave-equa- real. This is the von Neumann stability condition, the 
tion solutions, we will obtain the solutions of (3.1) directly is also the well-known restriction on the difference equé — * 
without elimination. tions for hyperbolic systems, which insures that the solt- M 
‘ 
bas 
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Fiaure 1.—Two lattices for solution of the primitive equations. 


tions of the difference equation and the solution of the 
differential equation coalesce in the limit At, Ar—0.* 


Taking this limit in (3.3) or (3.4) 
u,(t)=u,(0) cos v¢+h,(0) sin vit 
hilt) =h,(0) cos vt—u,(0) sin vit, 


(3.5a) 
(3.5b) 


the frequency equation being simply »;=ck;. 

Suppose now that we specify as initial conditions the 
values of u(x,0) at 2/Ar=0, 2, 4,..., 2p, and on the 
same line t=0 the values of h(z,0) at the intermediate 
points 2/fAz=1, 3, 5,...,2p—1. To initiate the lattice 
of figure la, we must move from h(z,0) to A(x,}At); if we 
do this by the uncentered step 


h(x,4At) =h (2,0) —408,u(2,0), 
the expansion coefficients h,(#At) will be 
hy(4At) =h,(0)—u,(0) sin 4v,At. 


(In these and in subsequent relations the bar is used to 


identify quantities associated with the lattice of figure 1b.) 
Substitution in (3.3) yields 


u,(t) =u,(0) cos vt+h,(0) (sec sin vt (3.6) 
hi(t) =h,(0) (sec 4v,At) cos vt —u,(0) sin (3.6b) 


Similarly, to initiate the lattice of figure 1b, we must 
rey from u(z, 0) to u(x, $At); proceeding as above, we 


“Convergence (coalescence) may be obtained for special choices of initial data even 
though ©>1 (see, for example, Platzman [5]). This may seem paradoxical because it 
Suggests that in the limit, the domain of dependence does not span the interval inter- 
cepted by two characteristics; but the special choices in question all involve analytic 
lultial data, capable of continuation. On the other hand, if one adopts the definition 
of convergence employed by Lax and Richtmyer [3], the condition «<1 is also necessary 
(a8 Well as sufficient), since they require convergence for arbitrary initial data. In their 
basic study, Lax and Richtmyer have established equivalence between stability and 
convergence, under very general conditions. 


=u,(0) (see 4v,At) cos vit+h,(0) sin (3.7) 


hi(t) =h,(0) cos vt —u,(0) (sec 4y,At) sin (3.7b) 
In (3.6) the initial data u,(0) are concordant, in the sense 
that they belong to the lattice on which the solution is 
represented; on the other hand, h,(0) in (3.6) are discord- 
ant data. Similarly, with respect to the lattice on which 
(3.7) represents a solution, the w,(0) are discordant and 
h,(0) are concordant data. 

Using the exact solution (3.5) for comparison, one may 
characterize (3.6) and (3.7) by saying that, if all initial 
data originate on the same line (t=0), the concordant 
data are propagated without distortion (of amplitude or 
phase) while the discordant data suffer an amplitude 
distortion sec }y,At. 

It is important to bear in mind that in the solution (3.6), 
on the lattice of figure la, the values of u,(t) are to be 
taken only for ¢/}4t=0, 2, 4, . . ., while the values of 
h,(t) are for t/}At=1, 3, 5, . . .; the converse is true of 
(3.7), which applies to the lattice of figure 1b. One can, 
of course, represent the two solutions in a single formula. 
This can be done with the aid of the functions 


F(t) =4[1+(—)"]+ —(—)"] sec 
G,(t) =4[1 see 


where n=t/}At; evidently, when n is even, F;=1 and 
G,=sec 4$v,At, and when n is odd, F,;=sec }v,At and 


G,=1. On this understanding, the solutions (3.6) and 
(3.7) are contained in the single representation 
u,(0)F,(t) cos vit+h,(0)G,(t) sin vt (3.88) 
h,(0) cos vit —u,(0)G,() sin vt, ——(3.8b) 
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in which (3.8a) and (3.8b) represent, respectively, w,(t) 
and h,(t) when t/}At is even, or u,(t) and h,(t) when 
t/}4t is odd. This solution exhibits a form which is 
characteristic of all double-lattice representations, namely, 
a “carrier” wave of period At, upon which the true solu- 
tion is modulated. The carrier wave is imposed by the 
functions F,(t) and G,(t). 

In connection with numerical computation, since 
truncation errors are determined primarily by the size of 
the mesh (Az, At), and since the two lattices of figure 1 
have the same mesh size, little is gained by computing on 
both lattices. Indeed, the double-lattice representation 
(3.8) evidently is more involved than (3.6) or (3.7), and 
tends to obscure the simpler features of the latter. 

A final comment will be helpful here, concerning a relax- 
ation of the end conditions w=0. Suppose that at the 
ends z=0 and z=L the canal is not closed, but that the 
flow there can be regulated independently of the motion 
in the canal; let w(0,t) and u(Z,t) be the end velocities thus 
imposed. The representation (2.2) may now be amended 
to read 


u(z,t) sin k,x+-u(0,t) +[u(L,t) —u(0,d)]2/L, 


while (3.2) is unaltered. Substitution in (3.1) leads to an 
inhomogeneous system for u,(t) and h,(t). The funda- 
mental modes are identical to those obtained above, but 
a particular solution must be added now to incorporate 
the forcing which comes from the ends. The integral 
form of the continuity equation is 


—o[u(L,t)—u(0,t)], 


obtained by summing (3.1b) over all z-points where h is 
defined. Clearly, this states merely that the change of 
mean level is determined by the net volume transport 
through the ends. 

In particular, if u(0,t) and u(Z,t) are constants (inde- 
pendent of t) we have the following particular solution 


u=u(0,t) +[u(L,t) 
h=—[u(L,t) —u(0,6)Jet/L 


(Since z and ¢ enter here linearly, this is a solution with 
zero truncation error.) @ The fact that A is a linear function 
of t obviously results from the constant net volume trans- 
port imposed at the ends. Mathematically, this may be 
interpreted as a resonant solution of the inhomogeneous 
system for uw, and h; coming from the mode corresponding 
to 1=0, the frequency of which is zero and coincident there- 
fore with the zero frequency inherent in the constant values 
of u(0,t) and w(Z,t).* 


4. THE VORTICITY EQUATION 


The linearized one-dimensional vorticity equation is 


*A thorough discussion of resonance in the application of end conditions to the finite- 
difference vorticity equation has been undertaken by Birchfield [1}. 


MONTHLY WEATHER REVIEW 


AUGUST 1958 
og 
+Bv=0 


v= dp/Odz, 


where 8 is the Rossby parameter and U is the zonal wind 
component, assumed uniform and constant. 
For the present purpose it suffices to take the case p=: 


ot or (4.1) 
Using central differences 
06,6, (4.2) 


where o=UAt/Az and 6,, 6, are the difference operators 
defined previously, as in (2.1). To aid in identifying the 
lattice structure of (4.2), and for comparison with the 
primitive equations discussed in the preceding section, we 
will denote vorticity values at the points z/}Ar=1, 3, 
5,. . ., 2p—1 by the symbol w and those at the points 
z/4dx=0, 2,4,. . ., 2p, by the symbol ¢. On reference 
to the lattices of figure 2, we find that equation (4.2) is 
equivalent to the system 


(4.38) 
(4.3b) 


6,¢=—ciw 


The structure of these equations is identical to that of the 
primitive equations (3.1), with ¢ playing the role of wu andw 
that of h. The significance of this equivalence is that (4.3) 
must admit of both forward and backward wave propa- 
gation (as in the primitive equations), whereas the parent 
differential equation (the vorticity equation) permits 
propagation in one direction only; in other words, (4.3) 
has the properties of a continuous hyperbolic system with 
two families of characteristics, while only one family is 
allowed by the vorticity equation. The intrusion of this 
virtual set of characteristics, and the associated truncation 
errors, are the penalties imposed by the use of central 
differences. 

The preceding remarks can be made more specific, 4 
follows. If one eliminates w from (4.3), the result is 


which is the wave equation (2.1). In other words, the 
¢-field on the lattices of figure 2 satisfies the finite-differ- 
ence wave equation. The corresponding differential 
equation may be obtained merely by differentiating (4.1) 
with respect to ¢, and eliminating the resulting mixed 
derivative by using (4.1) again; one finds 


ot? 


We now ask the question: how do the solutions of (4.5) 
compare with those of (4.1)? ma 

To examine this question, one must be very explicit 
about the boundary and initial conditions. For (4.5), ¥ 
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Ficure 2.—Two lattices for solution of the vorticity equation. 


impose ¢ and O¢/dt on the line t=0, and ¢=0 on the ends 
z=0 and z=L for all t; the latter condition corresponds to 
total reflection. For (4.1), we impose ¢ on the line t=0 
and ¢=0 on the end z=0; the latter condition presupposes 
U>0 (which we assume henceforth), so that z=0 is a 
point of inflow for the vorticity. The solution of (4.1) is 
then, with reference to figure 3, 


inside ApBoA; 
¢(z,t)=0 outside Ap 


Thus, the initial values of ¢ on AjBo, as well as the inflow 
values of ¢ on the line AyB,A, . . . are simply propagated 
without change along the characteristics z—Ut=constant, 
which are the only characteristics admitted by (4.1). 

To obtain an explicit and appropriate solution of (4.5), 
we will assume that d¢/Odt on the initial line only is specified 
in accordance with (4.1). This means that initially, the 
motion is organized strictly as a forward wave propaga- 
tion. A disturbance at any point on the initial line, such 
as D (fig. 3), therefore moves out along the forward 
characteristic z—Ut=constant, until the boundary at 
t= L is reached. There, total reflection occurs and the 
disturbance returns along the backward characteristic 
t+Ut=constant. It follows that in the triangle AyByCy 
the solutions of (4.1) and (4.5) coincide, since no point in 
this triangle is accessible by reflection; in particular, 
f (P )=¢(D). On the other hand, every point in the 
triangle C,ByA,, is accessible from two points on the initial 
line; for example, ¢(Q)=¢(D)—¢(EZ). In this triangle 


therefore, the solutions of (4.1) and (4.5) do not coincide, 
since (4.1) requires ¢{(Q)=¢(D). Evidently, (4.5) permits 
the Point Q to be affected by the reflected disturbance 
originating at Z on the initial line. Similarly, in the par- 


allelogram C,A,C,B,, we have ¢(R)=—s(D) from (4.4) 
and ¢(R)=0 from (4.1), while in the triangle B,C,A, we 
have ¢(S)=¢(E£) —¢(D) from (4.5) and ¢(S)=0 from (4.1). 

The preceding discussion shows that only in the triangle 
A,BoCy do the solutions of (4.5) and (4.1) coincide.* The 
finite-difference equations (4.3) will, of course, exhibit 
the properties of (4.5). We turn now to a consideration 
of these equations. 

In an earlier study of equation (4.2) the writer showed 
that an arbitrary imposition of ¢ on. the boundaries leads 
to a “neutral” prediction system; that is, a system whose 
modes .are neither damped nor amplified, apart from 
resonance (Platzman [6]). In particular, if one takes 
¢=0 at the ends, then the problem posed by (4.3): is 
identical to that of “‘tidal’’ motion in a closed canal, the 
solutions of which have been given in detail in the pre- 
ceding sections; hence, if in the latter solutions we replace 
u by ¢ and h by w, we obtain the solution of (4.3). From 
(2.2) and (3.2) we have 


a) sin ke (4.68) 
= ka; (4.6b) 


if all initial conditions are specified on the line t=0, we 
get solutions of the form (3.6a) and (3.6b), or (3.7a) and 
(3.7b), for the coefficients ¢,(t) and w,(t), respectively. 

It is pertinent to note that in the limit Az, At-0 the 
finite-difference solutions ¢ and w determined in the 
manner just described will coalesce with the continuous 
solutions of the equivalent primitive (or wave) equations, 

*In the particular case considered, with f=0 on the inflow boundary, we have coinci- 


dence of the two solutions also in the triangles AoCoBi, CiA;B2, and so on, where both 
solutions give f=0. 


1063 avevsT 1958 

I | 

oh lw ww wi | 
the Xx 
the 
(a) 
nts 
| 
3a) 
3b) 
the 
dw 
3) 
ant 
its 
3) 
th 
on 
ral 
as 
4) 
he 
al 
1) 
ad 
5) 
) 
| 


Be 


C, 


Figure 3.—Characteristic diagram for comparison of equations 
(4.1) and (4.3). 


but will not coalesce with the continuous solution of the 
vorticity equation, and in particular, the fields ¢ and w 
will not themselves coalesce. This anomalous behavior 
of the central-difference vorticity equation is readily 
understood from figure 3 and the related discussion of 
(4.5) given above. 

One way to circumvent the difficulties encountered in 
the central-difference vorticity equation is to use one-sided 
differences. In (4.1), we may take forward differences 
in ¢t and (if V>0) backward differences in z: 


where c=UAt/Az. This makes 
¢(2,t+ At)= (1 o)¢(z,t) + of Az,t). 


The solution of this equation is obtained most readily 
by induction. If ¢(0,t)=0, one finds 


where n=t/At, on the understanding that in the sum- 
mation, ¢(z,0)=0 where 

The use of one-sided differences can be extended quite 
naturally to the two-dimensional vorticity equation, and 
in this context has been reported by Bolin [2] as having 
been tested by the Stockholm group. For computational 
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stability one must make o<1; hence, if variations of 7 
are considered, one normally will have <1 over mos 
of the field, in which case by inspection of (4.7) it is cleg 
that the solution is highly damped.* 


5. THE STABILITY OF BOUNDARY CONDITIONS 


In the earlier study the writer [6] proved that insts. 
bilities are introduced by boundary conditions in which 
the vorticity at an outflow point is determined by “zero. 
or first-order extrapolation’? from the concurrent yor. 
ticities at neighboring interior points. The boundary 
conditions of ‘zero-order extrapolation” (case IT in the 
previous work) may be expressed 


4Az,t) ’ 


where the bar signifies elements of the lattice of figure 2b, 
Since the two lattices were not considered separately, 
these conditions evidently introduce a coupling between 
the lattices; the result is a weak but exponential instability 
in the fundamental modes. 

We will show now that if the outflow conditions just 
described are applied in such a way that the lattices are 
not coupled, then the instability is removed.t The 
appropriate revision is 


)=$(L—4z, t). 


(5.1) 


(5.2) 


Our proof (of stability) will consist of a construction of 
the fundamental modes. We begin with the expansions 


¢(z, t) = sin K,r (5.3a) 
w(x, t) = Silt cos (5.3b) 


which are similar to (4.6). Irrespective of the values 
assigned to x; these satisfy the basic equations (4.3) 
provided ¢,(¢) and w,(t) are suitable linear combinations 
of cos v¢ and sin v7, where »; is the frequency deter- 
mined by (2.4). Further, (5.3a) satisfies the inflow 
condition ¢(0,t)=0. In (5.3b), 2 is an absolute constant. 
We now effect the determination of x; by forcing (5.38) 
to comply with the outflow condition (5.2); this makes 


sin x,L=sin x,(L—Az), 
the roots of which are given by 

(l—}) Az) 

(l=1, 2,3, . . ., p—1). (5.4) 
There remains only the demonstration that, with « 
in (5.4), the representations in (5.3) are complete. This 


* In a thorough study of the finite-difference primitive and wave equations, Mihaljs® 
[4] has pointed out that when o=1 in (4.7), there are no truncation errors. 


For the conditions of “first-order extrapolation” (case III in the previous work), U 


coupling the lattices does not completely remove the instability, but eliminates its & 
ponential character. 
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hinges upon the following identities: 
‘a 


sin sin xvx(t)) 


if 141’) 


sin «,2(f) sin 0 if 


if 


COS K;X(w) COS (w) 
J 


where z(¢) and 2z(w) signify values of z at ¢-points and at 
wpoints, respectively, the z-sums being taken over internal 
points only. The first set of identities leads to the 


inversion 


=2(p— 4) COS K;,2(w), 


and the second set insures that when returned to (5.3), the 
latter expressions of ¢,(t) and w,(t) do, in fact, reproduce 
t(z,t) and w(z,t). 

For the determination of Q in (5.3b), we note that, in 
view of (5.4), cos x,(L—4Azr)=0; hence, from (5.3b), 


Q=w(L—}Az,t), 


which means that at the point adjacent to the outflow 
end, w is constant in this scheme. This completes the 
construction of the fundamental modes associated with 
the “uncoupled” outflow condition (5.2); evidently, these 
modes are neutral. 

An alternative proof of stability can be given in terms 
of the “‘semi-discrete” equation 


cy 2 

obtained from (4.4) by replacing the space derivative, but 
not the time derivative, by a difference quotient (see 
Platzman [6]). Let ¢(t) denote the column vector formed 


from the p—1 values of ¢(2,t) at internal points, and 
define the matrix 


(5.5) 


| 0 0} 
—1 2-1 
0 —!1 
Az 
0 
2 


Then, with the boundary condition (5.2) at z=L and 
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¢(0,t)=0 at z=0, the equation satisfied by ¢(t) is 
> 
2 
-(<) At. (5.6) 
This has fundamental modes exp (+7Ct), where 
C=(U/Az)?A. 


Clearly, the semi-discrete system is neutral if and only 
if the latent roots of A are real and non-negative. 

Now consider the fully-discrete system formed from 
(5.6): 


At. 
This has fundamental modes exp (+7Dt), where 


4 (sin (5.7) 


Evidently, the fully-discrete system is neutral only if the 
semi-discrete system is neutral. We now make the 
observation that the semi-discrete system (5.5) or (5.6) 
is precisely that which describes the transverse vibrations 
of a discretely-loaded, stretched string with one end fixed 
and one end perfectly free (with respect to transverse 
displacements). The mass of the string is assumed 
negligible, and the loading consists of p—1 equal loads 
each of which is concentrated at an interior {-point 
z/Ar=1, 2, 3,..., p—l. The end at z=0 is fixed, 
while at z=Z we have a free end which is not loaded; 
the condition (5.2) then expresses the requirement that 
the lateral component of tension must vanish at s=L. 
Since the modes for the loaded string are neutral (as is 
well known), it follows that the latent roots of C are reai; 
hence, from (5.7) we see that (for sufficiently small Af) 
the latent roots of D also are real. 

The preceding analogue permits some insight into the 
nature of the instability associated with the ‘coupled’ 
outflow condition. In the semi-discrete system, we have 
in general 


2 


but the “coupled” condition is 
=o(L— Az,t), 


so we get the following end condition for ¢: 
U 


For the stretched string, loaded in the manner previously 
described, with a free end at z=L, this condition may be 
interpreted as expressing the balance (at z=) between 
the lateral component of tension and a frictional resistance 
proportional to the displacement velocity of the end of 
the string along its support. Such a system clearly is 
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OF  O(w,x) 
dw O(f,¥), 
y) (6.2b) 


6 
> 


Fiaure 4.—The distribution of {- and w-vorticities in the z, y-plane. 


damped. ‘The modes of the corresponding fully-discrete 
system therefore are amplified. 


6. THE 2-DIMENSIONAL VORTICITY EQUATION 


In figure 4 is shown the appropriate spatial distribution 
of g- and w-vorticities for the two-dimensional vorticity 
equation 


ot 


oz Oy 


The w-vorticities are, as before, out of phase with the 
f-vorticities, by an interval 4At in time and 4As in space 
(note that figure 4 gives the distribution of these variables 
with respect to z and y, but not t). One may therefore 
identify two distinct lattices: one in which the ¢- and 
w-fields are defined respectively at t/}At=0, 2,4, . . . and 
1, 3, 5, . . . ; and the other in which the {- and w-fields 
are defined respectively at ¢/}At=1, 3, 5,... and 
0, 2, 4, .... In spite of the nonlinearity of (6.1), it 
is possible to formulate the difference equations in such a 
way that the two lattices remain permanently uncoupled, 
as will now be explained. 

Assuming the velocity components to be derivable from 
a stream function, equation (6.1) may be written as a 
pair of equations 


u (6.1) 


where x is the stream function for the w-vorticities, and 
¥ is the stream function for the ¢-vorticities: 

w=V?x; (6.3) 
The procedure for replacement of (6.2) by appropriate 
difference equations is rather apparent and need not be 
elaborated here. In connection with (6.3) however, two 
procedures may be mentioned. If the z,y-axes and the 
interval }As are used in the finite-difference Laplacian 
(to obtain the usual five-point formula), then concurrent 
values of ¢, w, y, and x are required, and in this way the 
two lattices would be coupled every time a relaxation js 
performed (that is, when (6.3) is solved). On the other 
hand, if axes inclined 45 degrees to the z,y-axes are used, 
together with an interval 4As72, then coupling of the 
lattices is avoided. 

Numerical solutions of the barotropic vorticity equs- 
tion have been obtained in the manner just described, in 
connection with an investigation of the prediction of 
hurricane movement, and will be reported elsewhere by 
Birchfield. 
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BAROTROPIC DIVERGENCE AND VERY LONG ATMOSPHERIC WAVES 
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ABSTRACT 


The problem of spurious retrogression of very long waves by the non-divergent barotropic forecasts is shown to 
be the same problem discussed extensively by Rossby, Yeh, and Bolin. This difficulty is due to the failure of the non- 
divergent model to allow properly for the mutual adjustment of wind and pressure fields. The equation of continuity 
for a homogeneous incompressible fluid with an upper free surface, proposed as a remedy by Rossby nearly 20 years 
ago, removes much of the difficulty. Further improvements are obtained by inclusion of a tropopause in the manner 
adopted by Bolin. The results of a series of 10 test forecasts are shown in verification of the function of the divergence 


in a barotropic model. 


1. INTRODUCTION 


The inability of a non-divergent barotropic forecast 
model to describe properly the behavior of the longest 
atmospheric waves has resulted in large systematic height 
errors of the prognostic charts issued by the Joint Numeri- 
cal Weather Prediction Unit. This difficulty was correctly 
diagnosed by Wolff [7], who also devised a method for 
automatic empirical correction to the forecasts. This 
consisted of an enforced stabilization of hemispheric wave 
numbers one, two, and three, such that atmospheric waves 
of this size were not allowed to change either position or 
intensity during the forecast. The remarkable success of 
Wolff's method was due to the fact that these very-large- 
scale disturbances are in fact quasi-stationary. 

However, it has been our ultimate objective to include 
in the forecast models proper descriptions of the physical 
processes which actually occur in the atmosphere, wher- 
ever possible. In looking for the appropriate physical 
mechanism responsible for the very long waves, one may 
note that their positions are relatively invariant, suggest- 
ing that geographically related factors are responsible for 
their formation and position. ‘Terrain-induced vertical 
motion and differential surface heating appear to be likely 
causes, 

The question of the behavior of these waves, once 
formed, can be considered apart from the question of their 
formation. Abundant discussion of this problem appears 
in the literature. It now appears that we have paid too 
little attention to the work of Rossby which bears on this 
question. 


2. THEORETICAL BACKGROUND 


In his classic paper of 1939, Rossby [4] pointed out that 
application of the vorticity equation to a non-divergent 
atmosphere did not properly take into account the adjust- 
Using a homogeneous 


ment of pressure and wind fields. 


incompressible atmosphere of thickness D, he used the 
equation of continuity 
dw, 1dD_ (1) 


dy D dt” 
where u and v are the speeds in the z and y directions. 
This led to the introduction of a divergence into the 
vorticity equation, such that 


dy_9 OD 
dt Dot 
where 7 is the absolute vorticity about a vertical axis, and 
where the geostrophic approximation is used to set 
oD, 
equation giving the phase velocity, C, of the divergent 
waves, which can be written in the form . 


0, (2) 


equal to zero. He then arrived at a frequency 


= (3) 


where C yp is the phase velocity of the non-divergent waves 
and N is the hemispheric wave number measured at about 
45° N. The coefficient 2.5 applies to a homogeneous 
atmosphere. The above equation and the corresponding 
one for the group velocity of this type of wave have been 
thoroughly discussed by Rossby [5] and by Yeh [8]. Ina 
footnote to his 1939 paper, Rossby [4] mentioned that an 
analysis of a two-layer model having an active lower layer 
and an inert upper layer indicated that the divergence 
introduced by these considerations might reach values as 
4 dD 


In his 1945 paper Rossby [5] again emphasized, very 
explicity, the necessity for allowing for the mutual adjust- 
ment of wind and pressure by means of the above type 
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of divergence term. Phillips’ [3] two-layer model can be 
specialized to describe Rossby’s two-layer atmosphere. 
In the case where the motion in the lower fluid is parallel 
to the contours of the interface (V-Vh=0), and the hori- 
zontal motion in the upper fluid is negligible, Phillips’ 
equations reduce to 


oz fn Oh 
oh , f fn Oh 
n 


In the above equations, h is the height of the interface, 
z can be thought of as the height of an isobaric surface in 
the lower fluid, V and 7» are the velocity and absolute 
vorticity of the lower fluid respectively, and f and g are 
the Coriolis parameter and acceleration of gravity. The 
factor ¢ is given by tbe ratio of the density of the upper 
fluid to that of the lower fluid. From equations (4) and 
(5) and from the lateral boundary conditions where 
02/0t=dh/dt—0 it is evident that 


08S OF 
Equation (6) is consistent with a continuity equation for 
the lower fluid, equivalent to 


1 oz 


where K=h/z. In using equation (6) as a barotropic fore- 
casting equation, it would be sufficient to use a mean 
value of h. 

Bolin [1] made some computations with a barotropic 
model having a tropopause. This model was the same as 
those described above. He used a factor of 4 for K(1—e), 
without giving a description of how he arrived at such a 
number. His results were widely misinterpreted at the 
time (and in particular at the JNWP Unit), since this 
feature was presented as a cure for excessive anticyclo- 
genesis. That problem has been shown by Shuman [6] 
to be related to the geostrophic approximation. Never- 
theless, Bolin’s results, while showing no improvement 
over the non-divergent forecasts in some respects, clearly 
resulted in greatly diminished height errors. It is evident 
now that Bolin was really controlling the very long waves, 
and that the reduced spread of excessive anticyclogenesis 
was merely a consequence of the smaller values of the 
group velocity obtained from his prognostic equation. 
At the JNWP Unit, we were aware of the problem of 
retrograde, very long waves only after we began to com- 
pute on a hemispheric grid and after Wolff’s analysis was 
completed. This was partly due to the fact that the cen- 
ters of height error from this source were very near the 
boundaries of our previously used grid and were confused 
with boundary errors. A suggestion by Dr. Norman 
Phillips led to the re-examination of the barotropic 
divergence. 


AUGUST 1958 
3. THE PROGNOSTIC EQUATION 


Equation (6) was reformulated in a finite-differeng 
form as follows: 


In this equation, y is the stream function obtained from 
the 500-mb. height through the balance equation [6] 7: 
and J are the finite-difference Laplacian and Jacobiay 
operators, d/m is the distance on the earth between grid 
points, and y is the same as 1/K (1—e). The evaluation 
of « seems to be a doubtful matter, since the atmospher 
is not composed of two homogeneous incompressible 
fluids. It was finally decided to use the atmosphere as an 
integrating machine in order to obtain a realistic value of g. 
Since we interpret h as the height of the tropopause, z as 
the height of the 500-mb. surface and consider that 


(9) 


we can observe the relative variation of h and 2z in the 
atmosphere. An inspection of tropopause and associated 
500-mb. charts reveals » to have an approximate value of 
4.* 

The second Jacobian in equation (8) is the standard 
large-scale mountain effect, where a, the ratio of surface to 
500-mb. wind speeds, is set at 0.2, p, is the pressure at the 
surface of the ground, and p is 1000 mb. 


4. RESULTS OF COMPUTATIONS 


Equation (8) was programmed for the IBM 704 com- 
puter for the 1977-point hemispheric grid used by the 
JNWP Unit. Computation time averages about 35 
seconds per time step. 

A series of 48-hour forecasts was made for varying values 


of u on one initial situation, that of 0000 amr, February J 


15, 1958. This case was chosen because of the very large 
amplitude of wave number one. 

The first forecast was made with a value of u=!, 
corresponding to a homogeneous atmosphere having 4 
free upper surface. The results of this are presented it 
figure 1. It is evident from this figure that the uope 
free surface model varies greatly from the non-diverget! 
model in its treatment of wave number one. The diffe- 
ence is in the direction required to reduce error. Figure ; 
shows the corresponding difference map for p=4. This 
particular 48-hour period was characterized by sever! 
of the most severe cyclogeneses observed all winter, and 


*If one follows the derivation of the equivalent barotropic model as given by Chame 
and Eliassen [2] but retains the term 7 ae at the lower boundary (w being the vertic# 


velocity in a pressure-time coordinate system), the resulting equation is the same / 
equation (8) with the important exception that the factor u is less than 1, being the ratio“ 
the surface to the mean wind speed. In the light of the results presented here, this form 
Jation of the equivalent barotropic model appears to underestimate the divergence neo 
vary to stabilize the longest atmospheric waves. 
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FievrE 1—(A) Stream function difference between 48-hour forecasts made by the barotropic model with »=1 and ~=0. (B) Stream 


dard function error of the 48-hour non-divergent barotropic forecast. 


ice to 


it the 
it should not be expected that any barotropic forecast 


model could eliminate all large errors in this situation. 
Fourier analyses were made by Cdr. Paul Wolff on 
the 48-hour stream function forecasts in order to discover 
com- ® the treatment of the various low-frequency components 
y the of the motion by the model. The results are presented in 
t 35% figure3. The large changes of displacement with changes 
in » for wave number one, compared with the minor 
alues changes for wave number five, are in good agreement 
ualy ® with the frequency equation (3). 
large Figure 4 shows the errors of the 500-mb. height fore- 
casts as a function of the coefficient uw. It is clear that 
u=1, ® the major improvement in the forecasts results from the 
ng 4 & free-surface approximation. It is also clear that little 
ed i HF difference is made when uw is varied between values of 
ope 3 and 8. 
igen A series of ten 48-hour forecasts was made in order 
iffe- BH to get a more representative idea of the level of accuracy 
ire ® of the forecasts. Comparisons were made with the JNWP 
Ths # Unit's operational barotropic forecast model (referred 
vert! HS to hereafter as “O” model), which includes the very-long 
an’ J wave stabilization described by Wolff [7]. This was 


thought to be a more relevant comparison than a com- ay LN 

aid with the unaltered non-divergent barotropic ZK RLY 

ms (“N” model), since the relative Jevels of accuracy 48 He Forecast /Tropopause_(-=4) — Nondivergent 


ato’ | Of the height forecasts of the O and N models were already 


po Presented by Wolff. Also, unpublished results by Wolff Ficure 2.—Stream function difference between 48-hour forecasts 
and Carstensen indicate a reduction of about 19 percent made by the barotropic model with «=4 and »=0. 
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At 
98 Hr DISPLACEMENT ih 
( 
60, Ne! 60, Nel 60, N=2 60, 60, N=3 60, N#5 the 
Southern Wave Northern Wove Southern Wave Northern Wave is 
40 40, 40. 40, 40; 4 
i 
20, 20 
-20! T ~20. -20. Est aft 
40, -40 -40 -40 ~4 of 
a ~60. ~60. ~60, ~60, th 
| ~80. ~60. ~60 res 
“100; “100 “100! “100 | 
~120. ~120, un 
? Ficure 3.—Displacement of Fourier components of the stream function in 48-hour forecasts vs. the values of coefficient uw. Forecasts for 
made from 0000 emt, February 15, 1958. a 
th 
st 


in the 48-hour wind errors of a series of winter forecasts Bt 
by the O model as compared with the N model. Table! & la 
500 presents the results of the ten comparative forecasts — lo 
: 4 (“D” refers to the divergent barotropic model). ev 
The average height errors for the ten forecasts shows & ot 
slight but significant improvement over the O model by 
the D model. A correspondingly large reduction of 
400. error is obtained by the D model as compared with N 
; model, considering the previously published results of 
Wolff. In comparison with the O model, the wind 
forecasts by the D model are not much improved. The 
improvement of the D model over the N model is, however, 


300. 


significant. 
= 5. CONCLUSIONS 
. . 
i 200 The success of forecasts made by a divergent barotropic 
| 7 model suggests strongly that large-scale divergences 1 
TaBLE 1.—Resulis of comparative forecasts made with the different 
models, D=divergent barotropic model, O=operational barotropi 
100- model, N=non-divergent barotropic model 
Initial map RMS height error (ft.) | RMS wind error (&t.) 
Date GMT| D N D N 
00 283 277 479| 27.1| 
141 170 |....----| 15.4 | 16.6 |------- 
Ficure 4.—Root-mean-square height error of the 48-hour forecast 00| 142 16.4] 16.6 |------- 
for different values of the coefficient Forecast made from 205 |........ 19.7 | 20.5 |----- 
0000 emt, February 15, 1958. 
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the atmosphere inhibit the motion of very long waves in 
the manner described by the model. Additional weight 
is given to the argument by the fact that Rossby clearly 
and explicitly emphasized the function of the barotropic 
divergence many years before numerical predictions could 
be made to test his conclusions. 

The inclusion of a tropopause in the barotropic model, 
after Bolin, gives further improvement. The estimation 
of the proper values of the coefficient y is a difficult matter, 
and is subject to some empiricism. It does not follow 
that a minimization of errors of the forecast heights with 
respect to the coefficient can lead to a correct estimation 
of its value, since it is possible that an overestimation 
might tend to compensate for correlated but physically 
unrelated errors. Investigations of mountain effects and 
surface heating with the aid of divergent forecast models 
may throw light on the formation of these very long waves. 

It has frequently been stated that numerical forecasts 
cannot be considered to be in competition with other 
forecast systems with respect to height values (at 500 
mb., for example). The above results, together with 
those previously published by Wolff, suggest that such 
statements, being based essentially on experience with 
the non-divergent forecast equation, and the consequently 
large height errors arising from misplacement of the very 
long waves, can no longer be considered valid. A re- 
evaluation of errors of height forecasts by numerical and 
other forecast methods appears to be desirable. 
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LOWEST TEMPERATURE IN CANADA 


ANDREW THOMSON 
Director, Meteorological Branch, Air Services, Toronto, Ont., Canada 
September 17, 1958 


Regarding the article “On the Lowest Temperature on 
Earth” appearing in the January 1958 issue of the Monthly 
Weather Review, I have a few remarks to make which 
might clarify the situation regarding the lowest tempera- 
ture ever recorded in Canada. 

Previous to February 1947, the record Canadian low 
temperature was —78.5° F., observed at Fort Good Hope, 
N. W. T. in December 1910, as indicated in the article. 
This absolute low temperature has been exceeded twice 
at Snag, Yukon, on February 2 and 3, 1947. The corrected 
minimum temperature on February 2, was —80.1° and 
on February 3, it was —81.4°. The value of —81° F., is 
taken as the absolute lowest temperature ever recorded 
in Canada. 

Since the article states that in Canada two official 
sources report conflicting data, it should be pointed out 
that the Monthly Weather Map, which is prepared about 
three weeks after the close of each month, must be 
considered as a preliminary or provisional report. As 
outlined in the article, the minimum originally reported 
from the station was —83° and was obtained by extrap- 
olation. When the thermometers were received at these 
Meteorological Headquarters, several months later, and 
calibration tests carried out, an official value of —81.4° F. 
was set. Since February 1947 there have been no official 
temperatures in Canada lower than —80°. 

I hope that this information will clarify the alleged 
discrepancy in the official Canadian meteorological 
publications. 


COMMENTS ON “THE ERROR IN NUMERICAL FORECASTS DUE 
TO RETROGRESSION OF ULTRA-LONG WAVES” 
PHILIP CLAPP 
Extended Forecast Section, U. S. Weather Bureau, Washington, D. C. 
September 22, 1958 


The paper by Paul M. Wolff (Monthly Weather Review, 
July 1958) discusses some practical results of one of a class 
of empirical numerical models which have been suggested 
in recent years. I would like to point out the interesting 
family resemblance between these models, because all of 
them contain some of the more important effects of tropo- 
spheric divergence resulting from forced flow over moun- 
tains, large-scale heat and cold sources, and other complex 
factors which account for the special behavior of long 
waves. 

I would also like to take this opportunity to discuss a 


more general formulation of this group of models, becauy 
I feel that more rapid progress may be made if researchex 
in this field realize that different individuals are ind. 
pendently attempting to develop the same fundament,| 
idea. This generalization is obtained by borrowing » 
assumption from the perturbation theory of hydrody. 
namics. It is assumed that the actual wind field in mid. 
troposphere consists of a perturbation-component super. 
imposed on a large-scale quasi-stationary flow. As in the 
perturbation method, it is assumed that this quai- 
stationary flow, which will be called hereafter the “basic 
current,” is given; and we seek to predict the evolution of 
the perturbation. Since, as will be clarified later, we ar 
not immediately concerned with mathematical difficulties, 
the basic flow need not be confined to a straight zonal 
current with little or no shear. It may have almost any 
form related to the general circulation, provided its rate 
of evolution is small compared to that of the perturbation. 
Another requirement, characteristic of the simple models 
to be described presently, is that it must be chosen soas 
to minimize the divergence of the perturbation. Variow 
possible forms of the basic current will be discussed later. 

The group of models which have been proposed up to 
the present may be formed through use of the simplest 
possible vorticity equation containing divergence: 

of 


(1) 


The symbols have their usual meaning and the equation 
applies to the actual wind field, which by definition is the 
sum of the basic current and the perturbation. 

Equation (1) may also be applied to the basic current, 
whose time rate-of-change is assumed to be small: 


(2) 


Here, the subscript B refers to quantities measured in the 
basic current. 

Use is now made of the assumption that the divergence 
of the perturbation is small. This enables us to replace 
the divergence of the actual flow by that of the basic flow, 
so that equations (1) and (2) may be combined, giving: 

of 


This is the prediction model proposed by the writer [1] 
where the basic current was taken to be the normal seé 
sonal circulation. In this particular adaptation, the 


(Continued on page 327) 
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g an 
ody. ABSTRACT 
mid- 
uper- Networks of data are simulated by interpolating height and wind from a hypothetically “correct’”’ analysis at a 
1 the uniform array of points. The interpolated values are added to random numbers whose statistics are typical of non- 
Uasi- systematic errors of observation, and are then regarded as genuine data. Such artificially constructed “data’’ for 
wae several networks of different densities are analyzed independently. The differences between these analyses and the 
. hypothetically “correct” analysis are taken to be representative of the initial analysis error under conditions of 
Dn of varying station density. 
> are Numerical forecasts computed from the different analyses are compared with the forecast made from the “‘correct”’ 
ties, initial data. Several such comparisons indicate that initial analysis errors do not grow to an important degree so 
onal long as the spacing between synoptic reports confines error fields to a scale smaller than that of the synoptic dis- 
any turbances. However, with data spacing comparable with that over existing regions of poor data coverage, initial 
. errors are amplified two-fold in a 48-hr. forecast interval. 
rate Further experiments, in which the data are analyzed objectively to eliminate inconsistencies of analysis, are 
carried out. 
dels 
0 
ious 1. INTRODUCTION data and analysis problem at length and has shown how 
= The problem of resolving the initial synoptic situation diBerenoes °4 barotzopio tendency computations can be 
noe. ‘ : ; : attributed to differences in analyses. Best [2] has carried 
in all its pertinent detail has often been mentioned as a 4 arg a 
lest 2 wet i. through several barotropic forecasts with similar conclu- 
contributing factor where critical weather developments te 
have challenged the forecaster and found his efforts this t y on 
inadequate. Although forecasters justify requests for q 
(1) : numerical forecast errors arising from initial analysis 
more data (at least among themselves) on the basis of that. i i , : . 
this argument, additional data are more often made function data density 
‘ion | available in connection with some expanded commercial — The experiment described below > actually an 
the | program such as a new air route where there is direct acjunct 6 8 mre theoretical treatment y Thompson 
operational need for current observsiions of weather con- [3]. 
casts on such routes are also contingent on adequate data conitea:” Oh h inemadinte Gill d 
coverage in other areas. This plea would perhaps be ts” later or 
2) weighed more heavily except for the subjectivity in the nal. 
the argument. No matter how elegant are the forecast P : 
parameters involved or how precisely observational data 2. DESIGN OF THE EXPERIMENT 
are employed in their evaluation, the final result is usually ‘ 
4 subjective combination which leaves the contribution of The experiment was designed to indicate differences in 
C the avowed critical facters eomewhat in doubt. forecasts (hereafter called errors) arising from differences 
m™ This argument could be greatly strengthened if a meas- 1 initial analyses of the same synoptic situations. Such 
we of change in forecast accuracy could be related ® definition of error seems entirely fair. The objective was 
quantitatively to a change in data coverage.’ The to measure quantitatively the reproduction of barotropic 
@ H woblem would then be one of economics with the in- ‘recess based on plentiful data by barotropic forecasts 
creased cost being weighed against the increased value of based on limited data. The similar task involving the real 
1, HE * more accurate forecast. Newton |1] has discussed this atmosphere or even more sophisticated models is certain 
“- Se to be more difficult. The barotropic forecasts were pro- 
he ‘Preliminary results of an investigation by Maj. E. O. Jess, USAF, have come tothe duced using an octagonal grid of 1977 points covering most 


attention of the writer. His approach is much the same as in the present paper although 
t measures of forecast error are employed. The two studies agree in a broad 


of the Northern Hemisphere (fig. 1). 
Four different analyses of two separate synoptic situa- 
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Figure 1.—Data coverage, shown by dots, for the octagonal grid of 1977 points currently used in the JNWP Unit (A) and successively 
reduced arrays (B, C, D,). 


tions were provided as input data for the forecasts in the 
following manner. The 500-mb. data from two routine op- 
erational analyses were first interpolated at the grid points 
indicated by large dots in figure 1A. It was then postu- 
lated that these grid-point data were the true heights of 


the 500-mb. surface at the points in question. The nex 
problem was to try to reproduce these true grid-pollt 
values with subjective analyses made under realistic ope 
ational conditions but with 500-mb. data coverage being 
provided in various reduced arrays. The analyst attempt 
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ingthe simulation had access to completeroutine operational 
surface analyses which contained reasonably complete 
coverage over most geographical regions including ocean 
gereas. The question was thus reduced to one of upper- 
air data coverage. Maps for 0300 emt for April 3 and 5, 
1957 were selected partly because data decks for points 
shown in figure 1A were readily available. Also, for a 
limited sampling, it is perhaps more meaningful to choose 
such a period having flow patterns of moderate intensity 
rather than either winter or summer extremes. 

The grid-point 500-mb. heights (considered to be uni- 
formly spaced perfect data) were modified by the applica- 
tion of random normally distributed errors having a mean 
value of 50 feet. Geostrophic winds measured at the grid 
points from the basic analysis were likewise amended by 
similarly applying a non-systematic 10-knot average wind 
speed adjustment. The realistic but uniformly distributed 
data thus obtained were next plotted on blank charts in 
three reduced arrays as shown in figures 1B, C, and D. 
In some of the illustrations these data arrays are referred 
to as Maximum, intermediate, minimum, and sub-mini- 
mum for reasons which will become obvious. The same 
experienced analyst proceeded to analyze one series, start- 
ing with the most sparse data array first. He was pro- 
vided with two preceding analyses having the same data 
array and was permitted use of normal differential analysis 
techniques. Once the coarse array had been analyzed, he 
was given the denser array and repeated the analysis 
process. 

Figure 2 presents the corresponding final analyses for 
the data arrays shown in figure 1 for the case of 0300 eur, 
April 3, 1957. As one might suspect, the patterns look 
strikingly similar, at least in a superficial way. Indeed 
one might say that the differences appear trivial if the 
chart is to serve as a basis for a subjective forecast. Sig- 
uificant differences do exist, however, and their repercus- 
sions in numerical forecasts are important. 


3. COMPUTATIONS AND RESULTS 


All competing forecasts were made using the operational 
barotropic forecast procedure being employed at the time 
of the experiment, except that special precautions were 
taken to minimize boundary error contamination. In 
particular an investigative code devised by ’Arnason [4] 
was employed to remove virtually all boundary inflow- 
outflow wind components from each set of input data. 
The technique involves replacing each boundary height 
by the mean of all boundary values and allows the 
discrepancy to be “faired” in with a weight of .4 at the 
first internal ring of points, a weight of .1 at the second 
"hg, and zero correction inward. Such an adjustment 
may, on occasion, do harm in the already questionable 
boundary region but, on the basis of other tests, it appears 
‘0 prevent large effects from penetrating meridionally far 
into the grid. In any case, rather elusive, but sometimes 
Mportant, minor inflow-outflow differences have thus 
been removed at a slight sacrifice in realism which does 
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not enter into the results since the error is defined as the 
difference between forecasts. 

Berggren [5] has recently shown the importance of 
subjective opinion among analysts in the interpretation of 
identical plotted data charts. Allowing one analyst to 
perform all the competing simulation analyses for each 
case was considered to be the best practical way of 
alleviating this difficulty. In this way any systematic 
habits or model concepts held by the analyst should not 
cause large random elements to enter into the results. 
The greatest chance for such a discontinuity to enter would 
be between the given operational analysis and the three 
simulated analyses but this difference enters more or less 
equally into all three resulting comparisons. Actually 
the analyst engaged in the second case was available for 
only a limited period and was unable to complete the 
analysis for the data array of figure 1B. The analyst for 
the first case completed the second case. One would 
expect this nonhomogeneous effect to be small since the 
data coverage dealt with was closest to the maximum 
thereby offering the least opportunity for varying inter- 
pretations. The results do not appear over-sensitive to 
this difficulty but differences between the two cases will 
be discussed in this light in connection with wind error 
results. The distribution of geostrophic wind errors was 
selected as the measure of forecast skill because Thomp- 
son’s conclusions involve the wind error and also because 
this measure is perhaps of most interest to those making 
direct use of 500-mb. prognostic charts. The several 
mean distances between observations expressed by the 
data arrays were selected for two reasons: (1) They 
approximate familiar arrays which presently exist over 
limited regions, and (2) error fields thus produced have 
different characteristic wavelengths. 

The coverage expressed in figure 1B closely approxi- 
mates that now in existence over most of Canada. Figure 
1C corresponds rather well to the present Atlantic coverage 
if one excludes reconnaissance data and also excludes 
a vast area south of 30° N. between Africa and the 
West Indies. Figure 1C approximates the spacing be- 
tween present subtropical Pacific islands and corresponds 
in general to Pacific coverage without reconnaissance. 
The results thus indicate error levels corresponding to 
uniform hemisphere-wide data distributions which exist 
at present over certain segments of the hemisphere. 
From the standpoint of scale comparisons, the mesh 
lengths in figure 1B and 1D permit error fields having 
half wavelengths of approximately 360 n. miles and 720 
n. miles, respectively. The half wavelengths of many 
synoptic disturbances having important local weather 
anomalies are shorter than the mesh length of figure 1D. 
There are likewise very few macro-scale disturbances 
(beyond their embryonic stages) which are characterized 
by half wavelengths less than the mesh size of figure 1B. 
The present study thus provides some empirical evidence 
bearing on the first three items listed by Thompson |3]: 
(1) Forecast range, (2) Initial wind error, and (3) Differ- 
ence in scale between disturbances and error fields. 
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forecasts from each input data deck were then pl 


duced by the current operational forecast 
initial field of stream potential. The 24-hr. ani 


This forecast procedure converts the input heigh! 
by using the balance equation (c. f. Shuman [6]) 
produces the forecast internally from t 


figure 1. 


April 3, 1957 based on data from the grid array with the corresponding letter it 
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were adjusted to a constant value as explained above. 
b. Hemispheric barotropic 24-hr. and 48-hr. height 


To summarize the procedure, the results which follow 
were obtained from the mentioned analysis procedure 


through the following computational steps: 
a. Boundary heights contained in the input data decks 


- Fiaure 2.—500-mb. analyses for Case 1, 0300 amr, 
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- Figure 3.—48-hour forecasts made using the analyses of figure 2 as initial data. 

48-hr. height output punch card decks are obtained points and a map of the corresponding wind errors 
pro- by inverting back through the balance equation. were produced for each such comparison. Further, 
ode ¢. For each of the cases, the initial, 24-hr., and 48-hr. the wind errors were sorted by 10-knot intervals to 
ghts height field decks involving the data mesh size of indicate a frequency distribution. 
ant figure 1A were compared with the corresponding Figure 3 displays the 48-hr. forecasts resulting from 
tins decks involving the data mesh sizes of figure 1B, 1C, the corresponding initial charts of figure 2. Important 


an and 1D. A map showing height differences at grid phase -differences-are to be noted particularly in the 
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(Left) Map of figure 2B compared with 2A. 
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(Right) Map of figure 2D compared with 2A. 


Fiaure 5.—Changes in the Case 1 height error fields at 48 hours (in hundreds of feet). 


Ficure 4.—Height error patterns (hundreds of feet) in Case 1 initial charts. 
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LE 1—Wind error ective analyses. 
(Percent of 1988 pot specified ) 
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Initial 24-hr. forecast 48-hr. forecast 

| Wind 
t nots) Max- | Max- | Max- | Max- | Max- | Max- | Max- | Max- | Max- 

Inter | Min | 8-min | Inter | Min | S-min | Inter | Min | S-min 

CASE 1 
a 66. 2 52.5 42.0 68.6 55.2 41.5 55.4 43.4 32.7 
10-20..--- 20.2 36.5 42.9 27.7 35.6 39.4 34.5 37.0 39.3 
-30..--- 4.1 8.5 11.3 3.4 7.8 14.8 7.5 13.4 17.6 
30-40.---- 4 1.8 3.0 -2 1.2 3.9 2.2 4.9 7.1 
0 
CASE 2 

ee 54.3 48.9 36.6 57.4 56.5 37.2 48.7 48.0 33.9 
10-20...-- 39.2 40.8 43.3 34. 5 35.7 38.7 39.4 40.4 35.1 
-30....- 5.9 8.3 14.2 6.1 6.2 16.2 8.9 9.5 17.9 
30-+40..... -6 1.5 4.3 1.3 1.3 5.5 1.7 1.8 7.0 
-4 1.3 -4 -2 14 -9 .3 3.0 
+2 +2 -6 0B 1.6 
-2 
$3 
trough feature in the region near California. Wind errors 


are even more striking. Figure 4 illustrates height error 
patterns for the comparison between initial charts having 
data coverages shown in figure 1A and 1B and similarly 
for coverages in figure 1A and 1D. In figure 4 the cellular 
pattern of errors for A—B is constrained to be small in 
scale as compared to similar quantities displayed for 
A-D. The changes in the error fields at 48 hours can be 
seen from the similar charts in figure 5. 

The principal wind error results are given in table 1. 
Errors of each category are in terms of percent of total 
area computed on the basis of 1785 internal points. 
Wind errors were computed from the height error charts 
by considering the gradients diagonally across the grid 
squares of the basic 1977-point grid. This length cor- 
responds to about 4° of latitude and approximates the 
portion of gradient normally used in hand measurements 
with a geostrophic wind scale. The results may be dis- 
played graphically in a variety of ways. With such a 
limited sample it is perhaps better to suppress some of the 
detail in the frequency distribution of errors. A simple 
two-category breakdown is presented in figure 6. Per- 
centage of area with wind errors over 20 knots is plotted 
as a function of forecast duration for the three comparisons 
involved. Even for rather intense winter regimes, a 20- 
knot wind error represents a large portion of the actual 
wind at 500 mb. From figure 6 it is first of all evident 
that the areas of large error increase rather uniformly as 
the data become more sparse with the sub-minimum 
(Pacific-type) coverage producing almost three times as 
much error as does the intermediate (Canadian-type) 
Coverage. In addition, it is of interest to note that the 


minimum error is reached at 24 hours on the lower two 
curves whereas it occurs at the outset for the most sparse 


According to Thompson [3], “. 


data array. .. if the 
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Fieure 6.—Wind errors versus forecast range, Case 1, subjective 
analysis. Max=grid array A in figure 1, Intermdt=grid array 
B, Min=grid array C, and Sub-Min=grid array D. 


scale of the initial error field were appreciably less than 
that of the true fluctuations, the error might actually 
decrease for a while.”” The evidence here supports his 
analysis with the lower error curves actually diminishing 
at 24 hours. Further, the behavior of the upper curve, as 
relates to reality, implies that there are limited regions in 
the central Pacific, where, on occasion, predictability is at 
best extremely marginal. Case 2 produced a similar 
picture but with the lower curve much closer to the middle 
curve. Since the analysis with data corresponding to 
Canadian coverage was carried out by a different analyst 
one might speculate that different habits of analysis might 
be appearing. A comparison of the error charts for the 
two lower curves suggests rather that the analysis with 
Canadian coverage was smoothed in excess of the error 
tolerance in the data thereby eliminating some of the 
real small-scale error prescribed by the network. An 
average of the two cases is presented in figure 7. 

The same evidence can be presented so as to emphasize 
more directly the question of data coverage. Figure 8 
uses the same error parameter in the vertical coordinate 
and uses spacing between observations as horizontal co- 
ordinate. This treatment permits an estimate of error 
magnitude for any uniform data array. Unfortunately 
actual observations are not uniformly spaced except to 
a rough degree over limited areas. In some actual arrays 
one sees that iso’ated reports are called upon to yield 
detail that is clearly impossible. If one counts the num- 
ber of observations in the outlined area of figure 9, ex- 
cluding ‘reconnaissance data, and weights all of. them 
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Fiaure 7.—Wind errors versus forecast range, average of Cases 1 
and 2, subjective analysis. 
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Fiaure 8.—Wind errors versus data density, Cases 1 and 2, sub- 
jective analysis. Same error parameter as in figure 7. 


evenly, the mean distance between reports is about 540 
n. miles. This is an extremely generous gesture since the 
analysis could obviously be improved by a more uniform 
redistribution. Even so, a uniform array of this dimen- 


sion already is in a rather intolerable range of error as 
shown by figure 8. With the addition of the reconnais- 
sance reports the again generous equivalent uniform array 
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Ficure 9.—Data coverage in eastern Pacific area. 


has a data spacing of about 300 n. miles. The array, of 
course, still is not uniform and some of the reports are 
not strictly synoptic. By weighting reports by the area 
represented and including “. .. North America, the 
Caribbean Sea, the North Atlantic and most of the 
Pacific . . .”, Thompson [3] arrived at an average dis- 
tance of 610 n. miles as an equivalent mesh length for 
much of the Northern Hemisphere. The implication here 
is that if a certain data density is considered adequate 
for the resolution of synoptic features of a given scale, 
then more closely spaced data are partially redundant. 
Viewing from the standpoint of resolving the broad-scale 
features of the flow, Newton [1] cited several instances 
where a slight redistribution would greatly increase the 
value of reporting stations. The present study certainly 
supports such suggestions. 


4. FURTHER TESTS USING OBJECTIVE ANALYSIS 


The previously mentioned inhomogeneity in the sub- 
jective analyses of the second case and the large amount 
of chart work involved led to an attempt to enlarge the 
sample by using objective analysis. Accordingly a code 
was produced which formed an artificial input data tape 
to be used by the objective analysis program presently in 
routine use by the JNWP Unit. The data tape thus 
prepared contained data for the same arrays and modi- 
fied them in the same manner as in the subjective analysis 
cases. One minor difference involved the application of 
random errors to the separate wind components thereby 
introducing variations in direction as well as speed. 

The subjective analyst relies on data at other levels 
and a qualitative image from past data as a starting 
point. The technique of objective analysis used here 
relies more on a highly satisfactory quantitative image 
from past data—namely a 12-hour prognostic chart. In 
the present study, provision for a representative first 
approximation to the analyses posed a problem. In 
keeping with the spirit of the first cases, such a first ap- 
proximation to the analysis should also reflect the analysis 
problems of the reduced data array. ‘To avoid extensive 
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Taste 2.—Wind error fr distribution: objective analyses. 
(Percent of 1221 points with specified error.) 


Initial 24-hr. forecast 48-hr. forecast 
error (knots) 
Wad Max- Max-S- Max- Max-8- Max- Max-8- 
Min min Min min Min min 
CASE 3 
56.6 42.8 58. 4 43.2 56.8 37.8 
33.7 39.1 33.4 38.2 34.8 38.3 
ee 7.4 13.7 6.9 13.2 7.0 15.6 
1.6 3.2 1.2 3.4 1.2 5.5 
.6 1.5 .2 2.0 
CASE 4 
60.0 46.8 65.2 44.4 58.8 41.4 
31.0 35.6 28.8 37.0 33.3 38.7 
6.9 12.0 5.0 12.4 6.3 13.7 
1.5 3.8 -9 4.9 1.6 4.5 


iteration by machine to accomplish this goal, an alternate 
simulation scheme was employed. The operationally 
available 24-hour forecast from the preceding day was 
substituted as a first approximation. This somewhat less 
perfect image of the current upper flow pattern is con- 
sidered to be comparable to an indirect derivation from 
concurrent lower-level information insofar as placement 
of major features is concerned. Such a forecast chart 
still would be expected to contain details of shear and 
curvature not resolvable with the most sparse data array 
of the experimental coverage network. Accordingly con- 
siderable smoothing was applied before substituting the 
forecast chart as the simulated first approximation for 
the analysis, 

The resulting wind error information analogous to that 
in table 1 is presented in table 2. Case 3 employed some- 
what less smoothing than did case 4. In general the 
main results of the subjective analysis experiments are 
duplicated. The comparisons involving the simulated 
intermediate data array (Canadian-type coverage) were 
not carried out. Figures 10 and 11 are to be compared 
with figure 6. Here one sees the same pattern. Initial 
analysis errors from the sub-minimum array (Pacific) are 
roughly twice as large as are the errors involving the 
minimum (Atlantic) array. Also the errors are reduced 
in both cases at 24 hours for the lower curve whereas they 
are not reduced in the upper curves. One important 
difference seems to be that the error growth rate is not 
as large for the upper curves in the objective analysis 
cases. Finally it should be pointed out the degree of 
smoothing was deliberately fixed in a separate trial case 
80 that the scale of the error pattern initially was pre- 
dominantly controlled by the mesh size of the reduced 
data arrays. To this extent the objective analyses were 
guided by the subjective tests. However this does not 
alter the result that, once the first approximation pro- 
cedure is fixed, we see the same pattern of behavior as in 
the subjective cases when the data array is changed. 
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Figure 10.—Wind error versus forecast range, Case 3, objective 


analysis. Compare with figure 5. 
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Figure 11.—Wind error versus forecast range, Case 4, objective 


analysis. Compare with figure 6. 


5. SUMMARY AND CONCLUDING REMARKS 


Two series of subjective analyses and the resulting 
barotropic forecasts indicate error behavior in accordance 


with that specified in Thompson’s investigation. 


Spe- 


cifically a measure of the forecast_wind error is shown to 
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be a function of forecast range, initial error, and scale. 


valuable consultations and to Mr. J. Murray Mitchel, 
Further, such errors are shown to represent a very real 


Jr. for his assistance in chart analysis 


practical problem in some areas for which forecasts are 
presently being issued. In such areas the future success 
of more sophisticated multi-parameter models may 
depend to a marked degree upon improvement in upper 
data coverage. Additional tests with objective analyses, 
which were somewhat dependent upon experience from 
the subjective analyses, suggest the same pattern of 
error behavior. Pending adoption of a baroclinic model 
for routine forecasts, additional tests of a similar char- 
acter for the baroclinic case are planned. It seems in- 
tuitively obvious that the three-dimensional analysis 
problem, which requires proper phasing between pressure 
and temperature fields, places even greater requirements 
for adequate data coverage. 
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Weather Note 


' WORLD RECORD LOW TEMPERATURE 


and 2000 em7, on August 9, 1958. The Sovietskaya minimum was 
measured with a remote electrical resistance thermometer with 4 
known correction, exposed in the shelter. ‘The temperature was 
lower immediately afterward but no calibration for thermometer,” 
the message stated. 


Antarctica has again broken the world low temperature record. 
A message received from Morton J. Rubin at Mirny, Antarctica 
announced the following minimum temperature readings: At Vostok 
(78°27 106°52’ E.) —85.8° C. (—122.4° F.) was measured be- 
tween 1200 emr August 7 and 0000 cmr August 8, 1958; at Soviet- 
skaya (78°24’ S., 87°35’ E.) —86.7° C. (—124.1° F.) between 1900 


CORRECTION 


Monraty WearuErR Review, vol. 86, July 1958, p. 253: In the second equation in column 
one should read 
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ESTIMATING THE VARIANCE OF DENSITY AT CONSTANT HEIGHT 
FROM CONSTANT PRESSURE DATA 


CONRAD P. MOOK 
Diamond Ordnance Fuze Laboratories, Department of the Army, Washington, D. C. 
(Manuscript received February 4, 1958; revised August 11, 1958] 


ABSTRACT 


A method is presented for utilizing atmospheric statistics developed for constant pressure surfaces in estimating 


the variance of density at constant height. 
lation of atmospheric statistics. 


1. INTRODUCTION 


The determination of the variation with altitude of 
the standard deviation of density is of interest in many 
problems in aeronautics and ballistics. Dines [1] was 
perhaps the first to present observational upper-air data 
showing the magnitude of the variability of density about 
the mean at certain fixed altitudes. His data extended 
upward to 13 km, 

However, in recent years, the practice of compiling 
upper-air statistics in terms of constant pressure surfaces 
rather than at intervals of constant altitude has become 
almost universal among the meteorological services of the 
world. This, though admirable from the point of view 
of synoptic meteorology, presents certain difficulties to 
those interested in obtaining summarized data for con- 
stant heights. The purpose of this report is to outline a 
procedure for approximating the vertical distribution of 
the standard deviation of density at constant height levels 
from constant pressure level statistics. 


2. THE VARIABILITY OF DENSITY 
AT FIXED ALTITUDES 


In 1919 Dines |1] showed how the standard deviation of 
the density at a given altitude could be estimated from 
the standard deviations of pressure and temperature. In 
Dines’ derivation, the equation of state is differentiated 
logarithmically to obtain: 


(1) 


Writing d, p, and t for the percentage variations of the 
quantities and squaring the result, he obtained: 


(pt) (2) 


which for small values of each term in equation (1), be- 
comes 


The method depends upon recent advances in the study of the interre- 


Therefore the problem of determining the standard 
deviation of density expressed in percent at constant alti- 
tude becomes that of determining the percentage stand- 
ard deviations of pressure and temperature at constant 
altitude and their coefficient of correlation. However 
as stated previously, data are now being made avail- 
able only in terms of constant pressure surfaces. Table 
1 contains data of this type for Goose Bay for January 
1952. 

The initial problem which we set out to solve in this 
study is to determine from the data shown in table 1 the 
corresponding standard deviations of pressure and tem- 
perature at constant altitude surfaces. A method of 
converting information, such as that shown in table 1, 
into the required form will be presented. The resulting 
standard deviations of density expressed in percent are 
then compared with similar statistics obtained by direct 
interpolations of density from constant pressure surfaces 
to fixed altitude levels. 


3. DETERMINING STANDARD DEVIATION OF PRESSURE 
AT CONSTANT ALTITUDE 


Buell [2] stated recently that “it seems intuitively 
obvious that the parameters used to describe the state of 
the atmosphere should be highly variable together or 
relatively constant together.” Klein [3] has prepared 
maps in which standard deviations of pressure at fixed 
levels are used interchangeably with standard deviations 
of height of constant pressure surfaces, after applying a 


TABLE 1.—Aver temperature and height for constant pressure 


surfaces and their standard iations, Goose Bay, Labrador, 
January 1952. 
P oT; 
(mb.) (° A.) (m.) 

256. 35 7.33 1319. 27 65. 37 
248. 60 6.05 | 3903. 53 105. 58 
226. 70 4.61 7648. 24 169. 50 
218. 82 5.99 | 11256. 57 164. 65 
BePaonqvcccypoysseegquocoresocoessseos 215. 59 4.45 15654. 14 14. 78 
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constant conversion factor based on the hydrostatic 
assumption. The factor to be used can be conveniently 
obtained from the Smithsonian Meteorological Tables 
(Table 60) [4]. 

Thus the height change statistics in table 1 may be 
converted to pressure change statistics using the factor 
appropriate to the mean temperature at the given pressure 
surface. For example, for an assumed mean temperature 
of —20° C. at 600 mb., the ratio is 8.1 mb. per 100 m., 
and therefore the standard deviation of pressure at con- 
stant level in the vicinity of 600 mb. becomes: 


p= 8.44 mb. (4) 


4. DETERMINING STANDARD DEVIATION OF TEMPER- 
ATURE AT CONSTANT ALTITUDE 


Stidd [5] has shown that the well known correlations 
between pressure and temperature at constant height 
(first presented by Dines [1]), also exist as a relationship 
between height and temperature at constant pressure. 
Therefore if a given temperature change in a constant 
pressure surface is associated with a change in height of 
that surface, the relationship can be used to determine 
approximately the temperature change which would occur 
at constant elevation.* 

The coefficient of correlation r,z¢ may be computed 
_ directly or may be inferred from o, and a7 statistics using 
the method of Stidd [5]. Both are shown in figure 1 which 
describes Goose Bay data for January 1952. 

The average departure from mean temperature asso- 
ciated with a departure of one meter in the height of the 
constant pressure surface is, at 600 mb.: 


or 


b=Prar 043 (5) 


where oy, is the standard deviation of temperature at 


constant pressure. 

Therefore at constant pressure, a 100-meter departure 
of height from the mean is associated with a 4.3° de- 
parture in temperature (of the same sign). Interpreted 
in terms of an increase in height of 100 meters, the in- 
crease in temperature which would occur at constant 
height would be approximately .65° C. greater than that 
at constant pressure, assuming a normal rate of decrease 
in temperature with altitude, as specified in the U. S. 
Standard Atmosphere in the troposphere [7]. 

Therefore 


4.95 
or (6) 


or 


where a7 is the standard deviation of temperature at con- 
stant altitude. 
At 600 mb. (from table 1) 


or=1.15 (6.05) =6.96 (7) 


*Since the submission of this paper, the author has seen a paper by Crossley [6] dis- 
inverse problem. . ‘ . 
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Fiaure 1.—Correlation coefficient between temperature change 
and height change of a constant pressure surface as computed 
and interpolated from and statistics. 


5. INTERPRETATION OF CORRELATION BETWEEN 
TEMPERATURE AND HEIGHT 


The remaining factor needed to evaluate the right side 
of equation (3) is rpr. A formula for computing the 
coefficient of correlation where z and y are departures 
from their respective mean values is 


(8) 


In the case wherein z and y are departures of height and 
temperature from their respective mean values, their 
standard deviations as well as their departures from the 
mean may be converted to constant level variations of 
pressure and temperature as shown in sections 3 and 4. 
Since the same multipliers are applied to both numerator 
and denominator of the right hand side of the formula for 
computing r, the two coefficients of correlation are equiv® 
lent, and thus may be used interchangeably. 


6. APPLICATION OF THE DINES EQUATION 
Having given the necessary variables in the Dines 
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squation we may proceed to the computation of the 
standard deviation of percentage density at Goose Bay 
at a level in the vicinity of 600 mb. as follows 


+ 4] 78)( Gas 6 


=1.95 percent 


This result may be compared with a computation of 
the standard deviation expressed as a percentage of mean 
density by actually interpolating to a fixed level between 
600 and 550 mb. Such a procedure yielded 1.8 percent. 

The procedure illustrated for the 600-mb. surface may 
be applied similarly at other elevations by using the 
following formula which combines the several steps ex- 


pressed in general notation: 


Here a is the factor obtained from the Smithsonian 
Meteorological Tables for converting o, to cp, T, is the 
assumed lapse rate (of NACA Standard Atmosphere, as 
used here), and other symbols are as previously defined. 

The results of computations for Goose Bay at other 
pressure surfaces for January 1952 are shown in table 2. 

As in the 600-mb. computation, a, the pressure change 
at constant level corresponding to a 10-meter height 
change of a constant pressure surface, was again obtained 
from the Smithsonian Meteorological Tables as shown in 
table 3. The temperatures assumed in table 3 are not 
precisely the mean values given in table 1, but they are 
within the limits defined by the standard deviation of 
temperature. This small departure from mean tempera- 
ture has little effect on the final result. For example, if 
at 200 mb. the assumed temperature is varied + 6° C., 
the resulting change in density variability is less than 0.1 
percent. 


TaBLe 2.—Standard deviation of percentage density for various 
_ Pressure surfaces at Goose Bay, Labrador, January 1952 data 


Pressure (mb.) Estimated Interpolated 3Xx100 
(%) (%) 
_ 23 19 
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tions. From [4] 


Assumed 
Pressure (mb.) tem a (mb./10 m.) 
—20 . 81 
.16 


In the computations, the results of which are shown in 
table 2, an isothermal lapse rate was assumed to exist in 
the vicinity of 200 mb. and 100 mb., assuming that these 
surfaces were in the stratosphere [7]. Under this assump- 
tion the standard deviations of temperature at constant 
level were further assumed to be equivalent to those at 
constant pressure. 


7. CONCLUSION 


A procedure has been developed for estimating the 
vertical distribution of the standard deviation of percent- 
age atmospheric density at constant height using data 
summarized with reference to constant pressure levels 
without specifying in advance the altitudes for which the 
estimates are to be made. However from these data, 
curves of percentage density variability vs. altitude may 
be constructed from which the percentage variability may 
be interpolated at any desired altitude. 
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THE WEATHER AND CIRCULATION OF AUGUST 1958’ 
A Month with an Unusual Temperature Reversal 
CHARLES M. WOFFINDEN 


Extended Forecast Section, U. S. Weather Bureau, Washington, D. C. 


1, INTRODUCTION 


Over North America the mean flow pattern during 
August 1958 (fig. 1) was characterized by pronounced 
blocking in the extreme north and by a planetary wave 
in relatively fast westerlies at middle latitudes concisting 
of a ridge in the western United States and a trough in the 
East. In contrast to July 1958 [3], these features of the 
westerly circulations were near their climatologically pre- 
ferred locations for the month and hence closely resembled 
the normal map but with slightly greater amplitude. They 
were responsible for an unusually warm regime in the 
West, which completely reversed the pattern of the pre- 
vious month. A July-August reversal of this magnitude 
has occurred only once before in the past 17 years. 


2. GENERAL CIRCULATION 


The blocking High over northern Canada and its deep 
companion depression across the pole accounted for the 
two largest departures from normal (+410 ft. and —500 
ft., respectively) which appear on the mean chart for the 
month (fig. 1). These two features became established 
early in the month, moved but little throughout the period, 
and completely dominated the weather of the polar basin 
and adjacent land areas. The polar Low was particularly 
persistent and vigorous, with a departure from normal 
at one stage in excess of 700 ft. (5-day mean map for Aug. 
16-20, fig. 5). However, although a pronounced depres- 
sion existed at high latitudes, the circulation over the 
remainder of Europe and Asia was weak and flat. 

The trough normally found in the western Pacific was 
weak and ill-defined, with the principal Pacific trough 
lying well to the east between 150° and 160° W. This 
corresponded roughly to its position during July [3] and 
illustrates the tendency of this regime to persist, even 
though the long-period climatological statistics indicate 
that troughs are not at all favored in this longitudinal 
band in midsummer [7]. Though the amplitude of this 
feature remained small, its center averaged 240 ft. below 
normal, and it maintained sufficient vigor to sustain the 
downstream ridge in western North America and effec- 
tively fill the northern portion of the west coast trough, 
leaving the latter effective only at lower latitudes and well 
offshore. This constitutes a significant change from July 
[3] when this depression, though similarly confined to 
southerly latitudes, was deeper and lay inshore along the 


1 8¢e Churts I-XVII following p. 328 for analyzed climatological data for the month. 


California coast. As a result of this filling and retro. 
gression, the cyclonically curved flow, which characterized 
the July pattern in the West, was replaced by anticyclonic 
streamlines in August which resulted in profound changes 
in temperature and precipitation anomalies. 

As mentioned, eastern North America and the Atlantic 
were strongly influenced by blocking. Twin warm Highs 
appear on the monthly mean pattern (fig. 1), one over 
northern Canada and the other over southern Greenland, 
These are associated with depressed polar vortices in 
eastern Canada and Iceland and the maintenance of a 
broad band of strong westerlies across the Atlantic at 
temperate latitudes. The jet stream, as illustrated by 
figure 2, reached a maximum speed of 14 m. p. s. at 
45° N., 50° W., or 4 m. p. s. higher than the normal at 
that point. This broad westerly wind stream and the 
deeper than normal trough along the east coast, combined 
to favor recurvature of tropical storms at sea followed by 
rapid motion across the Atlantic. 


3. HARMONIC ANALYSIS OF THE MEAN MONTHLY 
CHART 


In recent years increasing study has been given to the 
spectrum of amplitudes and energy of meridional motions 
on the scale of the general circulation. Graham [5], for 
example, studied the latitudinal variation of the 17,800-ft. 
contour of the 500-mb. surface and found that the sum of 
the first three harmonics adequately describes the normal 
map for January. He also found that the total of these 
components remained nearly stationary on daily 500-mb. 
charts during the month of January 1952. More recently, 
others [4, 6, 11] have dealt with energy spectra of the 
meridional component of the circulation. In general it 
has been found that two distinct wave number bands exist 
in the atmosphere, one incorporating waves one through 
four, and the other wave number five and above. The 
first class of waves is semi-permanent and may be ass0- 
ciated with the fixed topographic features of the earth's 
surface, while the second is connected with the smaller 
migratory systems. 

In order to study the amplitude spectrum of meridional 
flow for August 1958, a harmonic analysis of the 700-mb. 
flow pattern was prepared, following the procedutt 
described in [10]. The results are shown in figure 3 i 
which amplitudes are plotted against wave number fot 
the latitudes 75° N., 60° N., 45° N., and 30° N. The pre 
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Figure 1.—Mean 700-mb. height contours (solid) and departures from normal (dotted) both in tens of feet with troughs indicated by 


heavy lines, for August 1958. The circulation pattern over North America included a ridge in the West and trough in the East. 
The trough along the west coast was weak and effective only at low latitudes. 


eminence of the first harmonic at high latitudes is strik- 
ingly illustrated by the curve at 75° N. (and to a lesser 
extent by that for 60° N.). This reflects the existence of 
the two large anomaly centers on either side of the pole 
already described and indicates the degree of eccentricity 
introduced by such a circulation. The normal circulation 
for August similarly shows a maximum amplitude in wave 
one at 75° N. but only of about one-third the magnitude. 
At 60° N. the maximum amplitude was still associated 
with the first harmonic, with a weak secondary peak 
ning to appear in the fourth harmonic. This second 
maximum was also in evidence at 45° N. and 30° N. but 


displaced toward the higher wave numbers (dashed line). 
This characteristic dual maximum at temperate latitudes 
agrees roughly with the data of Eliassen [4] who found 
maximum amplitude at wave number one and again at 
wave number five, although his data were for a different 
level (500 mb.) and period (October 21 to November 30, 
1950). Conspicuous in the graph for 30° N. was the very 
large amplitude accompanying the low wave numbers. 
However, interestingly enough, the normal also exhibits a 
similar distribution, except that the amplitude of the first 
harmonic is even larger, a surprising 167 ft. The phase 
angle is such in each instance that it is largely attributable 
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Ficure 2.—Mean 700-mb. isotachs (in meters per second) for 
August 1958. Solid arrows indicate the position of the primary 
jet axis which roughly paralleled the northern border but dipped 
southward over the Atlantic because of blocking to the north. 


to the Asiatic monsoon and illustrates the marked eccen- 
tricity which this phenomenon impresses on the low- 
latitude circulation in midsummer. 

The first four harmonics accounted for 91 percent of the 
variability of the pattern for August 1958, and the addi- 
tion of the next four increased this figure to over 99 per- 
cent. For the normal chart for August, the first four 
harmonics similarly contain most of the variance (97 
percent), in agreement with the work of Graham [5], and 
only an additional three are needed to increase this to 99 
percent. Thus, waves one through four strongly dominate 
both the pattern for the month and its normal. This 
highlights the semi-permanence of these long planetary 
waves as opposed to the random nature of the shorter 
transient waves. However, the former are just those 
components which are inadequately handled by simple 
numerical models, which confirms the necessity of a differ- 
ent approach to this problem, as recently emphasized by 
Burger [1]. A step in this direction was taken by the 
Joint Numerical Weather Prediction Unit at Suitland in 
an attempt to take into consideration the special behavior 
of these ultra-long waves [2, 12]. Also it was this con- 
sideration which led Namias [9] to incorporate the extrac- 
tion of “latitudinal anomalies’ as an integral part of his 
“basic current” model. 


4. RETROGRESSION OF 5-DAY MEAN FEATURES 


As is often the case, considerable variation about the 
monthly mean occurred, and one of the more interesting 
trends discernible on 5-day mean charts was a long-period 
retrogression. To illustrate this tendency, three 5-day 
mean maps centered one week apart were selected (fig. 4). 

In thé initial state, represented by the chart: for August 
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Figure 3.—Amplitude in tens of feet of the zonal harmonic wave 
components at 75° N., 60° N., 45° N., and 30° N. for the 700-mb. 
mean chart for August 1958 plotted against wave number. Fo 
wave number one the amplitude was a maximum at both hig! 
and low latitudes, the former reflecting the eccentricity of the 
polar vortex and the latter the eccentricity imposed by tle 
Indian monsoon. A secondary peak appeared at 60° N. in the 
fourth harmonic which was displaced toward higher wave 


numbers. 


7 to 11 (fig. 4A), a fast westerly flow characterized the 
circulation with two pronounced troughs in the Pacifi 
and a third off the United States east coast. One weel 
later (fig. 4B) the pattern in the Pacific underwent a majo 
transformation. The deep trough off Japan filled ab- 
ruptly, resulting in an extended wavelength between th 
eastern Pacific trough and its new upstream counterpal! 
over eastern Asia. It adjusted to the new situation by 
retrograding to the position sbown in figure 4B, leaving” 
its wake a cut-off and greatly weakened trough near thé 


California coast. 
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America. This retrogression of the east coast trough 
coincided with the period when tropical storm Becky 
was approaching recurvature and permitted that storm to 
recurve off the coast in approximate conformity with the 
normal track for August. 

The Pacific trough continued its retrogression during 
the third week (fig. 4C), and once again the downstream 
features over North America followed suit. By this time 
the Canadian ridge was sufficiently far west to become 
realigned with the intensified maritime ridge in the 
eastern Pacific rather than with the High in the Plateau 
States as previously. Asa result, the United States trough 
shifted all the way back to the Mississippi Valley, which 
led to the only outbreak of really cold air into the central 
portion of the country. At this time, the Atlantic circula- 
tion also reacted by reversing the tendency for eastward 
motion of the Atlantic ridge. In fact, comparison of the 
change in position and strength of the Atlantic anti- 
cyclone from figure 4B to figure 4C indicates a net west- 
ward displacement of some 10°. This sharp retrogression 
of both the United States trough and the Atlantic ridge 
was instrumental in influencing hurricane Daisy to recurve 
farther west than either of its predecessors. 

It is noteworthy that, despite the overall retrogressive 
tendency described above, the trough in the eastern 
Atlantic succeeded in making eastward progress as far as 
the European coast (fig. 4C). However, the resultant 
overextended wavelength was subsequently adjusted by 
rapid discontinuous retrogression of this feature as well. 


5. TROPICAL STORMS IN RELATION TO 
THE MEAN CIRCULATION 


During August, three of the four storms which developed 
over tbe tropical Atlantic became hurricanes. Except 
for Ella, these fortunately remained at sea and recurved 
offshore, in rough accord with the mean circulation for the 
month (fig. 1). The first tropical storm, Becky, was 
detected near the Cape Verde Islands on August 10, 
though its exact position could not be established until a 
reconnaissance flight fixed the center about 20° east of 
San Juan on the 12th. From that point the storm con- 
tinued moving westward at a rapid rate, averaging over 20 
knots over its entire course prior to recurvature and at one 
time reaching a forward motion of 25 knots. As it came 
under the influence of the mean trough along the east coast 
it turned northward and recurved well offshore (Chart X). 
Becky never attained full maturity as a hurricane but 
remained fairly weak over most of its long trajectory. 

On the other hand, Cleo, the second storm, which 
followed Becky by only two days, was already a severe 
hurricane when first explored by reconnaissance on August 
14 when about 750 miles east of the Lesser Antilles. 
Although reports were insufficient to establish a definite 
track, it is probable that this storm had its origin in a 
tropical depression to the southeast of the Cape Verde 
Islands on the 11th. If this is the case, it progressed 
lather rapidly (about 21 m. p. h.) over the first portion of 


WS 


aN 


OR 


4 
C. AUG. 21-25, 1958 / 


Fiaure 4.—Five-day mean 700-mb. charts (in tens of feet) for (A) 
August 7-11, (B) August 14-18, and (C) August 21-25, 1958. 
The filling of the deep mean trough in the western Pacific resulted 
in retrogression of the west coast trough the following week (B) 
and of the United States trough and Atlantic ridge the third week 
(C). After the first period the west coast trough was weak and 
confined mainly to low latitudes. 


its journey, but deceleration must have begun near 45° W. 
prior to the reconnaissance penetration of the storm since 
thereafter it moved more slowly and shifted to a more 
north-northwesterly course (fig. 5). 

The key transformation, accompanying this shift in 
trajectory from that of Becky appears to have been the 
eastward motion of the Atlantic ridge and the splitting of 
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Ficure 5.—Five-day mean 700-mb. chart (in tens of feet) for 
August 16-20, 1958, with track of hurricane Cleo superimposed. 
The Atlantic ridge moved sufficiently far east that Cleo was able 
to penetrate the subtropical High. 


the subtropical high pressure cell (fig. 5). The retrogres- 
sion, described earlier as affecting the Pacific and the 
United States, had not yet begun to influence the Atlantic. 
As a result the ridge succeeded in making eastward prog- 
ress, though in the face of a lengthening wavelength. 
Subsequently (fig. 4C) this ridge did retrograde sharply, 
but it remained far enough east just long enough for Cleo 
to penetrate the subtropical high pressure belt. The 
evolution of this development is of some interest. On the 
5-day mean map for August 7-11 (fig. 4A) the western- 
most cell of the subtropical high pressure belt was the 
stronger. One week later, however (fig. 4B), the eastern 
center had grown at the expense of the western one, and 
two days later (fig. 5) it was definitely predominant. It is 
quite possible that the storm itself contributed to the 
anticyclogenesis and thus influenced its own steering 
pattern. 


6. THE TEMPERATURE REVERSAL 


The outstanding aspect of the temperature pattern for 
this month was the pronounced reversal from July which 
took place over most of the western portion of the country. 
Ordinarily, the reverse is true and the transition from July 
to August has usually been characterized by marked sta- 
bility of the temperature pattern. In fact, in a study based 
on the period 1942 through 1950 and later extended to 
1957, Namias [8] found persistence to be higher between 
this pair of months than for any other pair during the year. 
In his paper, persistence was expressed in terms of the per- 
centage of stations at which the temperature did not vary 
by more than one class (out of five). For the change from 
July to August this statistic turned out to be 80 percent, 
the maximum for the year. This August, however, only 59 
percent of the stations fell within this interval. Since the 
chance expectance of this event is also 59 percent, no 
significant persistence occurred. During the past 17 


5 MONTHLY MEAN TEMPERATURE 
CLASS DIFFERENCE BETWEEN 
JULY AND AUGUST 1958 


Ficure 6.—The number of classes the anomaly of temperature 
changed from July to August 1958. A marked reversal in 
temperature took place over a broad area in the West and North 
Central Plains. 


years, only August 1947 had a lower persistence tendency, 
with 51 percent. 

In August 1958, warming of two classes or more (cross 
hatched in fig. 6) encompassed an extensive area from 
the interior of California to Illinois, with maximum 
change (four classes) over Nevada and western Montana. 
The departure from normal of the monthly mean tem- 
perature rose 9.5° and 9.2° at Billings and Havre, Mont., 
respectively, and 9.1° at Winnemucca, Nev., to cite the 
three highest. These are unusually high figures and high- 
light the unusual nature of the change which occurred 
this month. 


The key changes in circulation accompanying this re- 
versal in pattern have already been described in section 2, 
namely, the strong ridging in the West which changed the 
flow from cyclonic in July to anticyclonic in August and 
relegated the west coast trough to a minor role offshore. 

Temperatures along the Pacific coast continued well 
above normal for the fourth successive month. Record 
warmth for August was reported from several locations in 
Oregon and from Seattle, Wash. Tatoosh Island experi- 
enced the warmest August in at least 56 years, and record 
or near record temperatures were the rule at most Cali- 
fornia coastal locations. This prolonged period of warmth 
along the Pacific coast has constituted one of the most 
distinctive characteristics of the temperature regime 80 
far this year and has been linked to abnormally warm 
ocean temperatures in previous articles of this series [3]. 
While a similar effect in all likelihood continued this 
month, the heat wave in the West can be adequately ex- 
plained in terms of the large-scale ridging in the West, 
above normal sunshine, and southerly components of flow. 

In general, temperatures in the rest of the country 
remained nearer normal, ranging slightly above normal 
over the Appalachians and eastward and slightly below 
normal over the States bordering the western Lakes and in 
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Departure of Average Temperature from Normal for the Week Ending 
— Midnight, Ls. t, August 10, 1958 


U_ S Weether Bureau 


~te Departure of Average Temperature from Normal for the Week Ending 
Midnight, 1s t, August 24, 1958 


Ficure 7.—Departure of average surface temperature from normal 
(° F.) for the weeks ending at midnight (A) August 10 and (B) 
August 24. The abnormal warmth which prevailed during 
the first period, especially over the Northern Plains, was followed 
two weeks later by an intense cold snap affecting central portions 
of the country. 


the Ohio and Lower Mississippi Valleys. Departures from 
normal were small because a warm regime early in the 
month was largely cancelled in the mean by a marked 
outbreak of polar air in the eastern two-thirds of the 
country toward month end. Figure 7A was selected to 
represent the warm regime. During this week tempera- 
tures were unseasonably high nearly everywhere but, 
relative to normal, were warmest in the northern Great 
Plains and Southwest where temperatures ranged up to 
6° above normal. Maximum temperatures in the northern 
Great Plains rose into the 100’s on the 8th and 9th, and 
Bismarck, N. Dak. recorded 107° F., the highest there 
since 1949. For these States and for those bordering the 
Lakes, where June and July had been unusually cool, this 
was the first real taste of summer weather. The circulation 
for this period (fig. 44) was characterized by broad west- 
erly flow along the northern border with mainly anti- 
cyclonic curvature. However, it will be recalled from 


section 4 that the stage was about to be set for a large-scale 

retrogression of planetary waves on 5-day mean maps. 

This was accomplished by the week of August 18-24, and 

during this period (fig. 4C) the principal mean trough 

became entrenched over the Mississippi Valley, bringing 
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Ficure 8.—Number of days in August 1958 with surface fronts 
of any type (within squares with sides approximately 500 miles). 
The band of maximum frontal activity extended from the Central 
Atlantic States west-northwestward. 


an abrupt termination to the hot weather in the Midwest. 
The strongest and most widespread outbreak of cold air 
of the month was initiated at this time and penetrated 
to the Gulf of Mexico. Figure 7B depicts the thermal 
distribution which ensued. Temperatures dropped rap- 
idly, and Madison, Wis. recorded 37° F. on the 25th, the 
lowest temperature reached there in August during 
the 80-year period of record. Near-freezing temperature 
occurred in the region of the western Lakes, and early 
season frost was reported in the higher elevations in the 
Black Hills of South Dakota. Thus, after only a brief 
period of summerlike weather in the Midwest, conditions 
more typical of fall were suddenly reintroduced. 


7. PRECIPITATION 


In general the precipitation for August, while not nearly 
as heavy as tor June and July, was nevertheless adequate 
over most areas. However, the July pattern of heavy 
rains and flooding through the Ohio Valley and Central 
Plains continued to dominate early August. Flood stages 
maintained in a number of streams in Illinois and Indiana; 
some of these streams had been in flood since about 
mid-June. A record 30.0-ft. stage was reached at New 
London, Mo. as a result of rains in excess of 6 in. in the 
upper reaches of the Salt River. This exceeded the pre- 
vious high reached in June 1928 by 1.2 ft. Thereafter, 
however, these excessive rains terminated as the ridge 
moved into the Midwest and brought warm drying winds 
to the area. 

During this period, when the subtropical ridge was at its 
strongest across the Atlantic and the Southern States 
(fig. 4A), onshore winds and easterly wave activity brought 
tropical rains to the Gulf States. In fact, this type of 
activity was quite prominent during August, and a number 
of 5-day mean maps (figs. 4, 5) exhibit pieces of weak 
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low-latitude troughs in the Gulf or along the Mexican 
border. These also served to furnish a supply of Gulf 
moisture so that frontal activity, which was fairly frequent 
in the East-Central States (fig. 8), succeeded in releasing 
substantial rainfall. In this manner a band of above 
normal precipitation for the month was produced, extend- 
ing from Texas northeastward to the Central Atlantic 
States (Chart III-B). One of the heaviest of these frontal 
rains transpired in central and eastern Oklahoma on the 
20th when 5 to 9 inches fell locally and caused considerable 
flood damage. 

Frequent shower activity was widespread in the South- 
west and, as a consequence, rainfall exceeded twice the 
normal over an extensive belt extending from California 
to southern Oregon (Chart III-B). This precipitation 
resulted from a moist tongue which curved anticyclonically 
around the western High and was sustained by Gulf 
moisture which was transported across Mexico by a broad 
current from the east-southeast. 

During this period Blue Canyon, Calif. reported the 
greatest number of thunderstorms ever to occur in August; 
Winnemucca, Nev. experienced the wettest August of 
record; and heavy showers and even a few flash floods 
occurred over the southern Sierras and desert regions 
of California. On the other hand, Grand Junction, Colo., 
which lay in the dry tongue on the opposite side of the 
mean ridge, (fig. 1), experienced tbe driest August in 
the station record which dates back to 1892 

In the Northwest, which remained hot and tinder-dry 
for the entire month, a serious range and forest fire hazard 
developed. Scores of fires were reported in the mountains 
of Oregon with several reaching significant proportions 
before being brought under control. Also, lightning 


started many fires in the Sequoia National Forest and 
other brush and forest covered regions of California. 
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VERTICAL WIND SHEARS NEAR THE CORE OF THE JET STREAM OVER 
THE NORTHEASTERN UNITED STATES, AUGUST 1-2, 1958 


ROBERT O. COLE AND LLOYD W. CHAMBERLAIN 
National Weather Analysis Center, U. S. Weather Bureau, Washington, D. C. 


1. INTRODUCTION 


One of the most valuable uses of knowledge of vertical 
wind shears is their application to planning optimum 
fight paths for jet aircraft. In forecasting winds for 
aircraft, their use as an incremental term has tended to 
simplify the task of interpolating winds between standard 
constant pressure levels. Furthermore, from the stand- 
point of flying safety and passenger comfort, the observa- 
tion of extreme vertical shears has been helpful in identi- 
fying regions of possible clear air turbulence. 

On August 1, 1958 there was evidence that strong 
vertical wind shears existed near the core of the jet stream 
over the northeastern United States, and within 24 hours 
the magnitude of these shears had diminished con- 
siderably. Wind shears computed directly from rawin 
reports (see table 1) indicated values as large as 21 
kt./1000 ft. Shears of this magnitude greatly exceed the 
value of 10 kt./1000 ft. which has been regarded as 
critical enough to divert a jet aircraft from its optimum 
flight path [1]. The purpose of this article is to analyze 
the conditions under which these extreme vertical wind 
shears occurred and to demonstrate their effective useful- 
ness by integrating them into a vertical wind shear chart. 


2. RELATIONSHIP OF SYNOPTIC FEATURES 
AND JET STREAM MODEL 


At 0000 amr, August 1, the 250-mb. chart (fig. 1A), 
which usually represents the standard level nearest the 
core of the jet stream, indicated a Low over eastern 
Canada with a flat short-wave trough extending south- 
westward from the Gulf of St. Lawrence to southeastern 
Michigan. The wind flow over the northeastern United 
States was generally from the west or northwest and 
contained two jet streams which tended to approach 


TaBLe 1.—Oceurrence of wind shears, August 1-2, 


Time M tude| Height th 

Location Date (emt) ‘| (kt. /1000 ft.) 
ldlewild, 1 0000 -12 41-43 2000 
Ouibou, Male. 1} 2000 
Payton, 1 1800 15 28-30 2000 
Dayton, 1 1800 -9 30-32 2000 
Dayton,’ Ohio. 2 1200 9 32-40 


each other over New England. One wind speed maximum 
was located near Maine with speeds of at least 120 kt., 
and the other was near Albany, N. Y. with speeds of more 
than 100 kt. 

Both jet stream positions seemed reasonably well 
located with respect to horizontal temperature gradient 
and contour spacing. Associated with the northern jet 
was a steep tropopause slope, while the southern jet 
appeared to be related to the surface frontal system. 
The surface front at this time was quasi-stationary with 
weak minor waves and extended westward from the 
Atlantic Ocean to a position south of Idlewild, on Long 
Island, N. Y., and into central Ohio. In view of the 
upper confluent pattern superimposed on the surface 
front, it seemed unlikely that any strong development was 
possible at this time. 

In conjunction with this synoptic pattern, Caribou, 
Maine reported a maximum vertical wind shear near 
the core of the jet of —10 kt./1000 ft. and Idlewild reported 
20 kt./1000 ft. The strong shear at Caribou occurred 
over a 5000-ft. layer located above the jet stream, to the 
rear of the trough, and in relation to a steep tropopause 
slope. When compared to the jet stream models described 
by either Endlich [2] or Riehl et al. [3], the magnitude of 
Caribou’s vertical wind shear seemed quite valid, espe- 
cially since it was apparently associated with strong 
horizontal wind shears to the left of the jet. 

While the northern jet stream was either over or just 
south of Caribou, the southern jet was at least 120 miles 
north of Idlewild. As would be expected, there was no 
evidence of strong horizontal temperature gradients south 
of the jet stream, yet Idlewild reported a positive shear of 
20 kt./1000 ft. over a 1000-ft. layer beneath the peak 
wind and a shear of —12 kt./1000 ft. over a 2000-ft. layer 
above the peak wind. The short-wave trough at this 
time was about 250 miles northwest of Idlewild and the 
slope of the tropopause in the area was negligible. 

According to the jet stream model, the largest vertical 
wind shears are usually found just above the core of the 
jet stream and on the left or cold side in conjunction with 
strong horizontal wind shears and temperature gradients. 
Strong vertical shears have also shown some tendency to 
be associated with strong peak winds, but the correlation 
was found by Dvoskin and Sissenwine [4] to be rather low. 
When compared to the jet stream model, the magnitude 
of Idlewild’s extreme vertical shears seemed highly 
questionable. 
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Figure 1.—250-mb. contour pattern including jet streams and isotach analysis (dashed). (A) 0000 amt, August 1; (B) 1200 cur, Av 
gust 1; (C) 0000 emt, August 2; and (D) 1200 emt, August 2, 1958. The eastward progression of the flat short-wave trough is ind- 
cated with its influence on the southern jet stream. The strong horizontal wind shears appeared in conjunction with the largest 


horizontal temperature gradient. 


Although the horizontal temperature gradient was 
rather weak, a computation was made to determine a 
reasonable vertical wind shear value for Idlewild, using 
the thermal wind equation 


oT* 
oz fT* os 


where V is the speed of the westerly wind, g is acceleration 
of gravity, f is Coriolis parameter, 7* is virtual tempera- 
ture, 2 is the vertical coordinate, and s is the horizontal 
coordinate, positive southward. For an average tempera- 
ture gradient south of the jet of 1° C. per 125 miles at the 
250-mb. level, the vertical wind shear at Idlewild was 
computed to be a little more than 1 kt./1000 ft. Actually, 
the computation showed no vertical wind shears for the 


temperature gradient with respect to Albany since the 
temperatures at Idlewild and Albany were identical. 
For the temperature gradient with respect to Washington, 
the shear at Idlewild was computed to be 2.5 kt./1000 ft 

From the large difference between the theoretical and 
observed vertical wind shears at Idlewild, it seemed likely 
that the strong reported shears were related to either 
instrumental error, reporting procedure, or microscale 
phenomena. In any case, the shears were regarded & 
having little or no effective value for prognostic purpose. 

At 1200 cmt, August 1, the short-wave trough at the 
250-mb. level (fig. 1B) was still quite flat but appeared to 
be located along a line through Maine, west of Albany, 
and in the vicinity of Pittsburgh. The northern Je 
stream was in about the same latitudinal position with 
respect to Caribou, but the wind maximum had moved 
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on eastward to Nova Scotia. The southern jet stream 
had moved southward to an east-west line through 
southern Pennsylvania with a wind maximum near Pitts- 
burgh exceeding 100 kt. The horizontal temperature 

jient to the north of the southern jet was estimated 
at 3° C./180 mi. 

Albany’s rawin at this time presented an analysis prob- 
lem in properly locating the jet since its wind at the 250- 
mb. level was within about 5 kt. of that at Pittsburgh. 
However, in view of the horizontal temperature gradient 
north of the jet and the apparent streamline relationship 
between Buffalo, Albany, and Portland (fig. 1B), Albany’s 
wind was regarded as unrepresentative. 

At the surface (fig. 2) the front was still weak with 
minor waves but had shifted a little southward to a 
position through Virginia and Kentucky. Except for 
Maine and portions of northern New York, there was a 
wide band of cloudiness over the area; and rain prevailed 
over much of Pennsylvania and nearby States. The posi- 
tion of the rainfall area had the conventional relation to the 
jet stream position and helped confirm the upper-air 
analysis. 

Coincident with the synoptic features for 1200 cmt, 
August 1, a vertical wind shear of —18 kt./1000 ft. over 
a 2000-ft. layer above the peak wind was reported by 
Caribou. While this value seemed rather excessive in 
relation to a reported peak wind of only 72 kt., there was 
some continuity with the past 12 hours for a large shear 
to exist. In relation to the tropopause, this large shear 
occurred in a zone of steep slope which was estimated at 
75 ft. per mile. With respect to the jet stream, it was 
difficult to determine from surrounding data whether 
the jet was north or south of Caribou; but in any event, 
the horizontal wind shear appeared rather weak on both 
sides of the jet and implied a smaller vertical wind shear 
than that reported. Computation of a vertical wind shear 
at Caribou by the thermal wind equation gave a vaiue of 
only .86 kt./1000 ft. In view of these features, Caribou’s 
shear of —18 kt./1000 ft. seemed much too large and 
had little or no useful value. ' 

Along the southern jet stream there were no vertical 

wind shear values exceeding the critical 10 kt./1000 ft. 
The largest shears computed directly from rawin data 
were 8 kt./1000 ft. at both Dayton and Pittsburgh which 
were near the base of the short-wave trough. Idlewild 
and Flint, Mich. had wind shear values of 3 kt./1000 ft. 
and 5 kt./1000 ft., respectively. 
_ Although these shear values did not conflict with the 
jet stream model to any great extent, it was rather sur- 
prising that Pittsburgh’s vertical shear was not larger. 
The core of the jet stream appeared to be located either 
over or just south of Pittsburgh, the peak wind was 
102 kt., and the 250-mb. temperature gradient over a 
130-mi. distance north of the jet was estimated at 1° C./25 
mit. Nevertheless, a computed vertical wind shear value 
of only 5.4 kt./1000 ft. was obtained from the thermal 
Wind equation. The slope of the tropopause over Pitts- 
burgh at this time was negligible. 
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Ficure 2.—Surface chart for 1200 emt, August 1, 1958, which was 
characteristic of the entire period. The surface front and the 
weather pattern appeared to be related to the southern jet. 
Area of cloudiness is hatched; precipitation areas are stippled. 


At 1800 emt, August 1, wind reports from Dayton 
indicated a positive vertical wind shear of 15 kt./1000 ft. 
over a 2000-ft. layer below the peak wind and —9 kt./1000 
ft. over a 2000-ft. layer above the peak wind. Since it was 
assumed that there had been little or no change in the 
synoptic pattern since 1200 ear, it seemed rather unusual 
for Dayton to have such strong vertical wind shears, 
expecially in conjunction with a peak wind of only 62 kt. 
Another suspicious factor was that the peak wind was 
reported at the 30,000-ft. level when only 6 hours previous 
it had been at 42,000 ft. From these facts with respect 
to continuity and the jet stream model, it appeared that 
these large vertical shears were unrepresentative and had 
little or no effective value for forecasting. 

By 0000 emt, August 2, the short-wave trough at the 
250-mb. level (fig. 1C) appeared to be near the east coast 
of the United States, but the contours remained cycloni- 
cally curved to the north of the southern jet stream from 
the Midwest to the Atlantic Ocean. The southern jet 
stream position was along a line north of Dayton, over 
Pittsburgh, and south of Idlewild, while the northern jet 
was about 100 mi. northeast of Caribou. The wind flow 
pattern was still west to northwest, but wind speeds at 
Albany, Caribou, Portland, and Nantucket had decreased 
to 50 kt. or less. Winds along the southern jet remained 
fairly strong with a maximum of at least 120 kt. over 
southern Pennsylvania. 

The tropopause chart at this time showed an unusual 
ridge across Pennsylvania separating low tropopause 
heights over Ohio from those in New England. The slope 
of the tropopause normal to the jet was greatest between 
Idlewild and Albany with a value estimated at 50 ft. 
per mile. Parallel to the jet stream, the tropopause slope 
between Dayton and Pittsburgh was about 34 ft. per mile. 
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Figure 3.—A composite time graph which shows the random fluctuation of several wind field parameters at stations near the core of the 
jet stream. Upper panel compares the height of tropopause with the height of the peak wind (dashed line). The center panel shows 
the magnitude of the peak wind; and the bottom panel compares the maximum positive vertical wind shears with the negative (dashed 
line) wind shears. All wind values were taken directly from reported rawin data. (A) Caribou, (B) Idlewild, (C) Flint. 


Concurrent with these synoptic features, Pittsburgh 
reported a vertical wind shear of —21 kt./1000 ft. over 
a 3000-ft. layer directly above the peak wind of 119 kt. 

$0 This strong shear was compatible with both the strength 

“s of the peak wind and the isotach pattern, which suggested 
that the horizontal wind shears were approaching a maxi- 

36 

30 


~ mum value. In figure 1C, the horizontal wind shears 
near Pittsburgh and to the south of the jet stream were 
130 130 SPEED IN KNOTS estimated to be a little less than 25 kt./deg. lat., and those 
110 SPEED IN KNOTS uo \—f to the north of the jet at about twice that value. Riehl’s 
90 90 \ [5] computations for maximum horizontal shear for 
70 yy 70 straight or slightly cyclonic flow have shown similar values. 
50 50 

30 30 


Although the occurrence of this strong shear at Pitts- 
burgh appeared quite reasonable with respect to the jet 
we * : stream model, the magnitude of the shear was regarded 
as excessively large. Kochanski [6] has computed vertical 

/1000FT.| shear of this magnitude by using horizontal temperature 
gradients of 10° C./deg. lat. In this case, the horizontal 
# temperature gradient normal to the jet was only 1° C./26 
i- mi., and the theoretical shear value computed from the 
thermal wind equation was only 5.4 kt./1000 ft. Although 
this computed value appeared too low at first, it seemed 

| in good agreement with the value of 6 kt./1000 ft. obtained 

$ by averaging Pittsburgh’s vertical wind shears over 4 

.| 10,000-ft. layer above and below the level of maximum 
Aug.1,1988 Aug.2,1958 ‘ug.1,1988 Aug.2,1958 wind. Therefore, in view of this analysis, the extreme 


Fiours 4.—Composite time graph of wind field parameters for (A) vertical wind shear at Pittsburgh at 0000 emr, August 2, 
Dayton and (B) Pittsburgh. was assumed to be the result of either some mi 
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phenomenon of short duration or of errors in wind meas- 
uring equipment or reporting procedure. 

At the same time, Pittsburgh had also reported a positive 
shear of 8 kt./1000 ft. over an 8000-ft. layer beneath the 

wind. This value seemed compatible with the shear 
of 5 kt./1000 ft. reported at Dayton and 9 kt./1000 ft. at 
Idlewild. All of these vertical wind shears were in general 
agreement with the jet stream model and acceptable for a 
useful analysis. As a matter of interest, a comparison 
between theoretical wind shears and the above shears 
determined directly from rawin data showed the latter 
to be about 3 times larger than those computed by the 
thermal wind equation. 

By 1200 emt, August 2, the short-wave trough at the 
250-mb. level (fig. 1D) was off the east coast of the United 
States. The southern jet stream was in about the same 
latitudinal position as 12 hours earlier, but wind speeds 
had diminished considerably. As might be expected, both 
positive and negative vertical wind shear values had also 
diminished. With the exception of Dayton, most wind 
reports along the jet stream indicated vertical shears of 
3 kt./1000 ft. or less. Dayton reported a positive shear 
of 9 kt./1000 ft. and a negative shear of 7 kt./1000 ft. 
over 8000-ft. and 10,000-ft. layers, respectively. Al- 
though these shears were neither extreme nor in direct 
conflict with the jet stream model, it seemed curious that 
they should have occurred with such strength near the 
end of a declining jet and in an area of rather weak 
horizontal temperature gradient. 


3. COMPATIBILITY OF DATA 


Since it has been established that wind data are subject 
to considerable variation near the core of the jet stream, 
figures 3 and 4 have been prepared to summarize the 
fluctuations during August 1-2 and to compare some of 
the wind field parameters. At 1200 emr on both August 
1 and 2, Caribou (fig. 3A, upper panel) presented an 
analysis problem when the height of the tropopause 
appeared to be out of phase with the height of the peak 
wind. Actually, subsequent investigation disclosed the 
presence of two tropopause points, and both met the 
currently adopted definition. Since the jet axis had heen 
analyzed within 150 miles of Caribou (fig. 1B), this out-of- 
phase relationship appeared to conflict with the jet model 
oo [3] but seemed compatible with Endlich’s model 
2 


Endlich has found similar out-of-phase relationships 
between the height of the tropopause and that of the peak 
wind and has classified them as type “A” tropopauses 
that are characteristic of the anticyclonic side of the jet. 
Type “A” is described as a high cold tropopause in con- 
Junction with a lower, weaker tropopause that is about 
10° C. warmer than the higher one. The average height 
of the lower tropopause was found to be 4000 ft. below 
the average height of the peak wind. While Caribou’s 
upper-air data exhibited some similarities to Endlich’s 
type “A” tropopause, it seemed quite likely that the 
analysis problem could best be solved by arbitrarily in- 
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creasing the wind speed adjacent to the lower tropopause 
within the limits of instrumental accuracy. The differ- 
ence in speed at the two tropopauses at Caribou was 
about 5 kt. at 1200 emt, August 1, and 17 kt. at 1200 
emt, August 2. 

As shown in figure 3B, the height of the peak wind at 
Idlewild appeared to fluctuate over a layer of 8000 ft. 
or more from 0600 emt, August 1 to 0600 emt, August 2. 
At the same time, the height of the tropopause showed 
only minor variations which tended to confirm Endlich’s 
[2] conclusion that the correlation between the height of 
the peak wind and the tropopause is not too significant 
for prognostic purposes. Also a comparison of the peak 
wind (center panel, fig. 3B) and maximum vertical wind 
shears (lower panel, fig. 3B) at Idlewild showed little or 
no relationship. 

At Flint (fig. 3C), the height of the tropopause and the 
height of the peak wind appeared almost ideally related. 
By contrast to this ideal relationship, the vertical wind 
shears at Flint, though small, were incompatible with the 
model. Reiter [7] and others have found that the magni- 
tude of negative shears is slightly larger than that of posi- 
tive ones. During August 1-2, the magnitudes of Flint’s 
vertical wind shears were the reverse of those found in the 
model. 

The relative strength of the positive and negative shears 
at Dayton (fig. 4A) fluctuated in random fashion during 
the period and showed little conformity to the jet stream 
model. Another incompatible feature in Dayton’s wind 
occurred at 1800 emt, August 1. At this time, the height 
of the peak wind was reported at 30,000 ft., which was more 
than 10,000 ft. lower than previous continuity would sug- 
gest. As was the case at Idlewild, there appeared to be 
little or no correlation between the speed of the peak wind 
at Dayton and the magnitude of the vertical wind shears. 

Of all the winds near the core of the jet stream during 
August 1-2, those at Pittsburgh (fig. 4B) seemed to fit the 
model best. The magnitude of the negative shears was 
larger than that of the positive shears in two out of three 
cases, and the largest shear tended to occur with the 
strongest wind. The heights of the tropopause exceeded 
the heights of the peak wind, but the difference between 
the two values averaged more than 10,000 ft. which seemed 
quite large when compared to Endlich’s and Riehl’s 
models. 

Most of the unrepresentative wind reports found near 
the core of the jet stream were attributed to the inability 
of wind-measuring equipment to measure accurately strong 
winds at high altitudes because of low elevation angles. 
According to U. 8. Signal Corps tests [4] on GMD-1 wind 
equipment, the root-mean-square wind error was found to 
be about 18 kt. at an altitude of 40,000 ft. and an elevation 
angle of 6 degrees. Since numerous oscillations of wind 
speed with height are often found in a complete rawin 
sounding, an error of 18 kt. at either the top or bottom of a 
1000-ft. layer could have introduced a fictitious vertical 
wind shear of 18 kt./1000 ft. If both the top and bottom 
were in error, a random relationship could have yielded a 
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Figure 5.—Smoothed vertical wind 
profiles along jet stream. ‘“X” 
\ reported wind speed, with direc- 
1 tion in 10’s of degrees at right edge 
of each section. Dot shows maxi- 
mum wind (90 percent of peak) at 
middle of maximum wind layer de- 
fined by dashed line connecting 
points of 80 percent peak wind. 
32 Large B = barotropic. DAY = 
7 Dayton, Ohio; PIT = Pittsburgh, 

Pa.; IDL = Idlewild, N. Y. 
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Figure 6.—Vertical wind shear charts which represent the 3-dimensional aspects of a smoothed wind field. Maximum winds and isotachs 
for (A) 1200 emt, August 1 and (C) 0000 emt, August 2 are within 10 percent of the reported peak winds and consistent with the 


250-mb. pattern. The height of maximum winds and the averaged vertical wind shears (dashed lines) for (B) 1200 emr, August 1 
and (D) 0000 emt, August 2 show the topography of the jet stream. Strongest vertical shear values lie north of the jet stream and 
their 12-hour displacement suggests continuous eastward movement. 


maximum fictitious shear of 36 kt./1000 ft. Questionable 
vertical wind shears similar to these values occurred at 
both Idlewild and Pittsburgh and were observed near the 
40,000-ft. level in conjunction with peak wind speeds of 
90 kt. and 119 kt., respectively. 


4. ANALYSIS OF VERTICAL WIND SHEAR CHARTS 


In the vertical wind shear charts prepared experimen- 
tally at NAWAC, the unrepresentative wind features near 
the core of the jet stream during August 1-2 appeared to 
be smoothed into a systematic and useful analysis. The 
procedure of constructing these charts at NAWAC is based 
on a technique proposed by Reiter [7], with the main differ- 
ence being that the vertical wind shears are averaged over 
4 10,000-ft. layer above and below the level of maximum 


wind. Then isolines are drawn directly for vertical shear 
values instead of depth of shear layer values as proposed 
by Reiter. In this manner, the magnitude of average 
vertical shear can then be used as an incremental term for 
rapidly estimating wind speeds for almost any altitude 
between the 500-mb. and 100-mb. levels by subtracting 
the shear term from the maximum wind value. 

Smoothed vertical wind profiles for several stations along 
the southern jet stream for the period August 1-2 are 
shown in figure 5 to illustrate the procedure of analysis, 
In general the profiles were drawn subjectively to conform 
with the assumption that winds near the jet stream in- 
crease and decrease linearly and have one peak velocity. 
In figure 5 the maximum wind is shown as the point equal 
to 90 percent of the peak wind, and_the depth of the layer 
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of maximum wind is estimated as the distance between 
points on the profile where the wind is equal to 80 percent 
of the peak wind. In this manner any wind speed within 
the layer of maximum wind is assumed to be within 10 
percent of the maximum wind, 

In figure 5A the smoothing procedure elevated the 
height of the maximum wind which was placed in the 
middle of the layer of maximum wind. Another item of 
interest is shown in figure 5 C, E, and G, in which the 
presence of an unreported wind was assumed by the 
straight-line technique. Both of these features obtained 
from the smoothed wind profile agreed with continuity 
and seemed quite possible since rawin records are known 
to contain numerous oscillations in speed which must be 
interpreted and smoothed before being reported. 

Ignoring the unrepresentative wind at 43,000 ft. made 
the wind profile at Pittsburgh for 0000 emt, August 2 
(fig. 5F) almost classic. The winds varied linearly with 
height and had one peak velocity. The vertical wind 
shears appeared to be greater above the peak wind than 
those below. Vertical wind shears computed directly 
from rawin data were 9 kt./1000 ft. below and —21 kt./ 
1000 ft. above the peak. When smoothed and averaged 
over a 10,000-ft. layer above and below the level of maxi- 
mum wind, the magnitude of the vertical wind shears was 
reduced to only 6 kt./1000 ft. which was in good agree- 
ment with the 5.4 kt./1000 ft. shear computed from the 
thermal wind equation as previously noted. 

The wind profiles at Pittsburgh and Idlewild (fig. 5 G 
and K) for 1200 emt, August 2, showed a peak wind of 
less than 50 kt. Profiles which fell into this category or 
had a layer of maximum wind greater than 15,000 ft. 
were classified as ‘“‘barotropic” (B) since there was little 
or no vertical wind shear. 

Although the extreme vertical wind shears observed 
during August 1-2 had been rendered imperceptible by 
the smoothing process, the resulting values presented an 
organized three-dimensional picture of the wind field near 
the core of the jet stream when plotted on the vertical 
wind shear chart (fig. 6). The jet stream position for 
1200 emt, August 1 (fig. 6A) was in essentially the same 
position as that on the 250-mb. chart (fig. 1B). The 
isotach pattern near the core of the jet had been altered 
from that at the 250-mb. level, but then the winds had 
been computed with a 10 percent tolerance and were 
representative of a layer about 7000 ft. deep (fig. 5 B, 
E, I). 

The height of the jet stream core (fig. 6B) appeared to 
undulate from 40,000 ft. at Dayton to 37,000 ft. at Pitts- 
burgh and back again to 39,000 ft. at Idlewild. In 
relation to the synoptic situation, the lower heights of 
the jet stream core seemed to be associated with the posi- 
tion of the flat short-wave trough at the 250-mb. level 
(fig. 1B). 

Although there were no extreme vertical shears reported 
along the southern jet stream at this time, both Dayton 
and Pittsburgh reported modest and reasonable values of 
8 kt./1000 ft. These values were reduced to 4 kt./1000 
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ft. during the smoothing process. In relation to the syp. 
rounding data (fig. 6B) shear values of at least 4 kt./1009 
ft. appeared to exist not only along the jet axis in Ohio 
and Pennsylvania but also about 2° of latitude north of 
that position and suggested the possibility of at least » 
6 kt./1000 ft. shear in an area without data. This pattern 
was in substantial agreement with the jet stream modd, 

By 0000 cmt, August 2, the jet stream (fig. 6C) was in 
about the same position as that 12 hours earlier and was 
in good agreement with the 250-mb. chart (fig. 1C). The 
height of the jet core (fig. 6D) varied between 38,000 and 
39,000 ft., and followed continuity, with the lower height 
values appearing to the rear of the flat short-wave trough 
at the 250-mb. level. Strong horizontal wind shears were 
observed to the north of the jet maximum, and it was at 
this time that Pittsburgh reported an extreme vertical 
shear of 21 kt./1000 ft. 

During the smoothing process, this extreme shear at 
Pittsburgh was reduced to a value of 6 kt./1000 ft.; how- 
ever, there was now 12-hour continuity for a shear of the 
lesser magnitude. While this interpolated value of 6 kt,/ 
1000 ft. had moved eastward into Pennsylvania, the 4 
kt./1000 ft. isoline had progressed to Idlewild. Although 
the period of continuity had been relatively short, the 
extreme vertical wind shears near the core of the jet 
stream during August 1-2 seemed to fit a more usable 
pattern when smoothed than when allowed to fluctuate in 
random fashion. 

As previously discussed, the procedure of averaging wind 
speeds and vertical shears appeared to be justified on the 


TABLE 2.—Vertical wind shear (kt./1000 ft.) 


Greatest | Computed From 
>) from smoothed 
Station compu thermal Eat 
from rawin wind (20,000-ft. 
equation layer) 
0000 cmt, August 1 
1200 cmt, August 1 
Ji Jol OF Jo; -3 23 
-18 86 12 
0000 cut, August 2 
9 24 
1200 GMT, August 2 
Idlewild -2 3.5 
Pittsburgh. . 2 3.2 
Flint... 1 
Caribou. -3 23 
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basis that wind equipment in current use does not have 
sufficient accuracy reliably to detect small-scale phe- 
nomena in the upper troposphere. The microscale feature 
of clear air turbulence has been found [8] by numerous 
jet aircraft penetrations to have an average size of less 
than 50 mi. in length and less than 2000 ft. in depth. 
In. view of these considerations, it is believed more de- 
sirable to locate zones of turbulence with reference to a 
model rather than to unrepresentative wind data of 
questionable value. 


5. CONCLUSIONS 


Most of the extreme vertical wind shears computed 
from rawin data near the core of the jet stream during 
August 1-2 showed little or no agreement with the jet 
stream models, lacked continuity, and appeared unrepre- 
sentative when compared to surrounding stations. When 
the wind data were smoothed to remove the unrepresenta~ 
tive features and the vertical shears were averaged over 
a 20,000-ft. layer, a pattern emerged on the vertical 
wind shear chart that indicated some prognostic value 
in planning wind forecasts for jet aircraft. A comparison 
of the greatest vertical wind shears computed directly 
from rawin data with those computed from the thermal 
wind equation and those obtained from the smoothed 
wind profile is made in table 2. The largest average 
vertical wind shears were compatible with the jet stream 
model and were located north of the jet in conjunction 


with the strongest horizontal wind shears and temperature 
gradient, Although some errors were undoubtedly in- 
troduced in the smoothing process, they were considered 
insignificant in the broadscale sense. 
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Correspondence 
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divergence of the actual flow is replaced by that of the 
normal, which presumably is in turn a reflection of the 
e planetary influences mentioned earlier. 

Variants of this method, based on recognition that the 
second term on the right of equation (3) is related to the 
errors of the barotropic model, have been suggested by 
Berson and Williams [3]. 

The relation to Wolff’s model may be seen if the basic 
current is defined to be the sum of the first three har- 
monics of the actual flow pattern. The analogy becomes 
more exact if the fictitious changes in the three harmonics 
are removed at the end of each time-step. Perhaps some 
difference still remains due to the fact that Wolff’s model 
permits the changes (or tendency) of the three harmonics 
to be influenced by non-linear interaction with the shorter 
waves, while this is not permitted when using equation (3). 

A method very similar to that of Wolff, but dealing 


with one-dimensional wave motions, is that proposed by 
Graham [4]. 

Another in this family is the present operational model 
of the Extended Forecast Section, U. S. Weather Bureau, 
proposed by Namias [5]. Here a “fictitious” wind (Vp) 
is defined by the equation: 


Vr=V—(Va—V3) (4) 


where V is the actual wind, V, is the basic current at the 
same point and V¢ is the latitudinal average, or zonal 
component of the basic current for the latitude of V. 
The predictive equation is then obtained by assuming 
that the absolute vorticity of the fictitious current is 


(5) 
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The close correspondence between equations (5) and 
(3) can be seen by expanding them using the definition of 
the perturbation current (V*): 


V=V*+Vs (6) 


In this manner equations (3) and (5) become (7) and 
(8) respectively: 


(tet f)—V (7) 


It will be seen that these two equations are identical 
except that in equation (8) the zonal component of the 
basic current replaces its value at a point. Physically, 
the difference appears to be related to the fact that in the 
Extended Forecast Section’s operational model the per- 
turbation is initially uncoupled from the basic current. 
This appears to prevent any important non-linear exchange 
of energy between these two currents, while such an 
exchange is permitted using equation (3). 

The normal circulation was first used as the basic cur- 
rent in equation (4), but later it was decided to replace 
this with a large-scale flow pattern more characteristic of 
the particular season during which the forecast is prepared. 
A monthly-mean chart, centered on forecast day, is now 
used for this purpose. Tests indicate that this results in 
considerable improvement (Namias, to be published). 
Although the monthly-mean circulation evolves more 
rapidly than the normal, this apparently is more than 
compensated by a reduction in magnitude of the perturba- 
tion divergence. 

It must be clear that perturbation models such as those 
described above need not be derived solely by using the 
approximate vorticity equation (1). For example, the 
writer is at present working with one of several single- 
parameter models where a direct estimate is made of the 
divergence associated with large-scale planetary waves. 
In this case, equation (1) cannot be used in the perturba- 
tion method, and at least one more term must be added, 
preferably one which expresses some unknown or difficult- 
to-measure property of the actual flow. This unknown 
property is then replaced by the corresponding known 
property of the basic current. 
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It is believed that this procedure has the advantage that 
the models are constrained to behave in a fashion consist. 
ent with the observed behavior of the general circulation, 
because in essence the errors of any model are replaced by 
estimates based on experience. It is also felt that this 
procedure may lead to a better physical understanding of 
the circulation. In this connection it may be noted that 
the success of Wolff’s model, which makes use of the ob- 
served stationary character of long waves, stimulated the 
successful search for a new operational model which con- 
tains a more satisfactory physical basis for this charac. 
teristic [6]. This new model is based on the theoretical] 
work of Rossby [7]. Another important theoretical paper 
which treats the special problem of planetary waves is 
that of Burger [8]. Perhaps a more general formulation 
of the empirical approach, together with these theoretical 
studies, can lead to a practical solution of the preferred 
geographical locations and observed motions of the long 
waves. 
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Chart I. A. Average Temperature (°F.) at Surface, August 1958. 
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A. Based on reports from over 900 Weather Bureau and cooperative stations. The monthly average is half the sum of the monthly 
average maximum and monthly average minimum, which are the average of the daily maxima and daily minima, respectively. 
B. Departures from normal are based on the 30-yr. normals (1921-50) for Weather Bureau stations and on means of 
25 years or more (mostly 1931-55) for cooperative stations. 
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B. Departure of Average Temperature from Normal (°F.), August 1958. 
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Chart III. A. Departure of Precipitation from Normal (Inches), August 1958. 
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B. Percentage of Normal Precipitation, August 1958. 


Normal! monthly precipitation amounts are computed from the records for 1921-50 for Weather Bureau stations and from 
records of 25 years or more (mostly 1931-55) for cooperative stations. 
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Chart VI. 


A. Percentag 


e of Sky Cover Between Sunrise and Sunset, August 1958. 


B. Percentage of Normal Sky Cover Between Sunrise and Sunset, August 1958. 
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"A. In addition to cloudiness, sky cover includes obscuration of the sky by fog, smoke, snow, etc. Chart based on 
visual observations made hourly at Weather Bureau stations and averaged over the month. B. Computations 
of normal amount of sky cover are made for stations having at least 10 years of record. 
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Chart VII. A. Percentage of Possible Sunshine, August 1958. 
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B. Percentage of Normal Sunshine, August 1958. 
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A. Computed from total number of hours of observed sunshine in relation to total number of possible hours of 
sunshine during month. B. Normals are computed for stations having at least 10 years of record. 
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